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Chapitre 1 

Introduction aux équations de Lagrange 


1.1 Equations de Lagrange pour une particule 

1.1.1 Equations de Lagrange 

Considérons le cas particulier d’une particule astreinte à se déplacer, sans frottement, sur 
une courbe plane contenue dans le plan xOy. La courbe sur laquelle est astreinte à se déplacer 
la particule de masse m, est le lieu des points dont les coordonnées vérifient les relations : 

f z = 0 

\ f(x,y ) = 0 

La première relation correspond au plan xOy . La seconde relation représente l’équation de la 
trajectoire dans ce plan. Ces deux relations définissent les équations des liaisons appelées souvent 
liaisons. Le nombre de degrés de liberté est égal au nombre de coordonnées qui représentent la 
position de m (trois dans le cas général) moins le nombre de liaisons (deux dans le cas particulier 
étudié ici). La particule possède donc un degré de liberté. Il faut choisir une variable q pour 
repérer sa position. Cette variable est appelée coordonnée généralisée. Il est possible d’exprimer 
le vecteur position f de la particule en fonction de la coordonnée généralisée q par la relation : 
f = r(q). 

Soit F la résultante de toutes les forces agissant sur la particule. La relation fondamentale 
de la dynamique s’écrit : 

d 2 f dv 

F = rri——r = m— 
dt 2 dt 

où « = ^ est la vitesse de la particule. 

Soit 5W le travail fourni par la force F lors d’un déplacement infinitésimal ôf : 

SW = F -ôf 

Le déplacement infinitésimal ôf peut s’écrire en fonction de la variation ôq de la coordonnée 
généralisée q : 

df 

ôf= —ôq 
oq 

Dans ce cas le travail SW peut se mettre la forme : 


df 

ÔW = F ■ —ôq 
oq 
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Introduction aux équations de Lagrange 


On appelle force généralisée conjuguée de q, ou (/-composante de la force, la quantité F q 
définie par : 

„ SW -, dr 

q ôq dq 

Par conséquent SW s’écrit : 

SW = F q Sq 

En tenant compte de la relation fondamentale de la dynamique, cette expression peut égale- 
ment s’écrire : 


D’autre part : 


Sachant que 


on obtient 


dv dr 

ovv = m— ■ - — oq 
dt oq 


d 

j 

dr 

dv 

dr 

1 dl 

d 

dr 

dt 

V • 

dq 

dt 

'dq +V ' 

dt 

dq 


d 

dr 

d 

dr 

dt 

dq 

dq 

dt 


dv 

dq 


dv dr 
dt dq 


d _ dr 
dt 1 dq 



Le vecteur vitesse v, peut aussi s’écrire : 


dr drdq dr . 

dt dq dt dq ^ 


D’où la relation : 


dr dv 
dq dq 


dv dr 
dt dq 


d_ 

dt 




Sachant que 


d_ 

dq 



d_ 

dq 


2 V ' V 



et que 


d_ 

dq 



d_ 

dq 


1 _ _ 
-v ■ v 



on obtient 


dv dr 
dt dq 

L’expression du travail SW peut alors s’écrire : 


d 

' d 

1 2 

d 

1 2 

dt 

dq 

2 V 

dq 

2 V 


SW = m 


{- 

' d 

1 2 

\ dt 

dq 

2° J 


d_ 

dq 



Sq 
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Si on note T = \mv 2 l’énergie cinétique de la masse m , on obtient finalement : 


w - < 5 


dT 

dq 


-SI» 


On obtient finalement les deux expressions équivalentes du travail 5W 

' dT 1 


a i_ 
dt 


dq 


ST' 


On en déduit l’équation de d’Alembert pour un système à un degré de liberté : 


d_ 

dt 


dT 

dq 


_dT =p 
dq q 


1 . 1.2 Cas des systèmes conservatifs 

Dans les systèmes conservatifs, la force appliquée au système dérive d’un potentiel U et elle 
s’écrit : 

F -Æ 

q ~ dq 


L’équation de Lagrange devient alors : 


a l 
dt 


dT 

dq 


dT _ _dU 
dq dq 


dU 


Généralement l’énergie potentielle U ne dépend pas de la vitesse, c’est-dire que — — = 0. 

dq 


L’équation de Lagrange peut alors s’écrire : 


c l_ 
dt 


d(T-U) 

dq 


d(T — U) 
dq 


= 0 


On introduit la fonction de Lagrange ( ou lagrangien du système ) qui est la différence de 
l’énergie cinétique et de l’énergie potentielle : 

C = T-U 


D’où la forme de l’équation de Lagrange dans le cas d’un système conservatif : 


d_ 

dt 


dC 

dq 


dC n 
-~dq~ ] 


1.1.3 Cas des forces de frottement dépendant de la vitesse 
Equation de Lagrange 

Considérons une situation physique dans laquelle la particule est soumise à des forces de 
frottement de viscosité dont la résultante / est de la forme : 

/ = -av 


Pour calculer la force généralisée f q correspondante, nous utilisons la définition du paragraphe 
précédent : 


fq 


f- 


df 



dr 2 dq 
dq dt 
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Introduction aux équations de Lagrange 


Cette dernière expression peut se mettre sous la forme : 


fq = ~PQ 


avec 


P = a 


Ôr 2 
dq 


Si en plus des forces qui dérivent d’un potentiel il existe des forces de frottement de viscosité, 
l’équation de Lagrange s’écrit : 


d_ ïdT 
dt dq 


dT 

dq 


Fu,q + fq 


OÙ F Utq 


dU 

dq 


représente les forces qui dérivent d’un potentiel. D’où : 


j d_ r dC 

dt dq 


dC 

dq 


-Pq 


Fonction dissipation 

Calculons le travail SWf fourni par la force de frottement pendant un intervalle de temps St 
pour un déplacement ôr : 


SWf = f ■ ôr = —a v 2 ôt 

La quantité de chaleur 6Q gagnée par le système en interaction avec la particule, est telle 
que : 


ôQ = a v 2 ôt 


S(D 

Soit P d = -f- la puissance dissipée par les forces de frottement sous forme de chaleur : 
dt 

P d = a v 2 

Cette puissance dissipée peut être exprimée en fonction de q, par : 

I 2 


Pd = a 


dr 

2 

drdq 

dt 

= Oi 

dq dt 


= Pq 


■2 


Par définition, la fonction dissipation est égale à la demi-puissance dissipée : 

1 _ 1 


V=-P d =~Pq 


■2 


La (/-composante f q de la force de frottement peut alors s’écrire : 

dV 


fq ~ dq 


L’équation de Lagrange s’écrit alors : 

A l— 

dt dq 


dC dV_ 

~ dq + dq ^ 
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1.1.4 Cas d’une force extérieure dépendant du temps 

Considérons le cas plus général d’une force extérieure dépendant du temps agissant sur un 
système qui est le siège de forces de frottement qui dérivent d’une fonction dissipation V. Soit 
F eq la (/-composante de la force extérieure. Dans ce cas l’équation de Lagrange peut s’écrire sous 
l’une des deux formes équivalentes suivantes : 


d 

'dC 

de 

dt 

8q 

d q 

d ' 

dC 

dC 

dt 

d q 

Ôq 


= F eq - 13 q 


dV 

Ôq 


e,q 


1.2 Système à plusieurs degrés de liberté 


Dans le cas général d’un système à plusieurs degrés de liberté, il y a autant d’équations 
de Lagrange que de degrés de liberté. Ainsi, si le système possède N degrés de liberté, il est 
nécessaire d’avoir N coordonnées généralisées q % (i = 1, 2, ...., N) ; nous aurons ainsi N équations 
de Lagrange : 

d_ 
dt 


dC dC dV „ 

Wi “% + %“ eÆ (v " ’ } 


La qj— composante de la force généralisée extérieure est définie par : 


jp — 

r e,qi ~ 


SW 

Sqi 


Sqi^O 


Dans cette expression ÔW représente le travail des forces extérieures résultant d’une variation 
Sqi de la coordonnée qi telle que les coordonnées q 3 ^i soient constantes {Sqj^i = 0). 


1.3 Exercices 


Exercice 1 : On considère un point matériel astreint à se déplacer sur un cercle de rayon R et 
de centre O contenu dans le plan xOy. 

1. Traduire la liaison par une ou des relations mathématiques; quel est le nombre de degrés 
de liberté de ce point ? 

2. Quelles sont les coordonnées généralisées que l’on peut utiliser pour repérer ce point ? 

Exercice 2 : On considère un point matériel astreint à se déplacer sur une sphère. Répondre 
aux mêmes questions que l’exercice précédent. 

Exercice 3 : Pour repérer la position d’un solide dans l’espace, il faut repérer la position de 
trois points non alignés A, B et C de ce solide. 

1. Traduire les liaisons physiques par des relations mathématiques ; quel est le nombre de 
degrés de liberté de ce solide ? 

2. Quelles sont les coordonnées généralisées les plus couramment utilisées pour décrire le 
mouvement d’un solide ? 

3. Quel est le nombre de degrés de liberté pour un solide qui possède : 

(a) un point üxe? 

(b) deux points fixes ? 
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Exercice 4 : On considère une haltère constituée de deux masses identiques m, supposées 

ponctuelles, reliées par une tige de longueur a, de diamètre et de masse négligeables. 

1. Comment s’écrit mathématiquement la liaison entre les deux masses? 

2. Quel est le nombre de degrés de liberté de ce S 3 ^stème? 

Exercice 5 : On considère une masse M qui glisse sans frottement selon une droite sur un 
plan horizontal. Elle est reliée à un bâti fixe par un ressort parfait de raideur k, colinéaire avec 
la trajectoire. 

1. Quel est le nombre de degrés de liberté? 

2. Quelles sont les forces qui s’exercent sur la masse M. Quelles sont celles qui dérivent d’un 
potentiel ? Quelles sont celles qui ne travaillent pas ? 

3. Calculer l’énergie cinétique et l’énergie potentielle de ce système ; en déduire l’équation 
différentielle du mouvement par la méthode des équations de Lagrange. 

4. Etablir l’équation différentielle du mouvement en utilisant la seconde loi de Newton ; que 
remarque-t-on ? Quelles sont les forces qui n’interviennent pas dans l’équation de Lagrange 
et qui sont prises en compte dans les équations de Newton ? Quelle est leur particularité ? 

Exercice 6 : On considère un pendule simple constitué d’une masse m reliée à un point fixe 
O par un fil de longueur i et de masse négligeable. Cette masse peut osciller librement dans le 
plan vertical xOy. 

1. Quel est le nombre de degrés de liberté de ce système? Quelles sont les coordonnées géné- 
ralisées les plus pratiques à utiliser ? Ecrire les coordonnées x et y de la masse m dans le 
repère xOy en fonction des coordonnées généralisées choisies. 

2. Quelles sont les forces qui s’exercent sur la masse m. Quelles sont celles qui dérivent d’un 
potentiel ? Quelles sont celles dont le travail n’est pas nul au cours du mouvement ? 

3. Etablir les équations du mouvement par la méthode des équations de Lagrange. 

4. Ecrire les équations du mouvement par la méthode de Newton ; retrouve-t-on le même 
résultat que par la méthode de Lagrange ? Déterminer le module de l’action du fil sur la 
masse m ; pouvait-on déterminer ce module par la méthode de Lagrange ? Commenter le 
résultat. 

Exercice 7 : Etudier le mouvement d’un cylindre de masse M et de rayon R, qui roule sans 
glisser le long de la ligne de plus grande pente d’un plan incliné qui fait un angle ip avec l’hori- 
zontale. 

Exercice 8 : Etudier à l’aide des équations de Lagrange, le mouvement d’une masse M qui 
glisse sur un plan incliné faisant un angle <p avec l’horizontale, avec un coefficient de frottement 
de glissement y,. La masse est soumise de plus à une force F(t) parallèle au plan incliné. 

Exercice 9 : Etudier, à l’aide des équations de Lagrange, le mouvement d’un cylindre de masse 
M et de rayon R autour de son axe de révolution fixé horizontalement, entraîné en rotation par 
l’action de forces extérieures dont le moment par rapport à l’axe de rotation est M{t). 

Exercice 10 : Une particule de masse m est lâchée sans vitesse initiale dans un fluide caractérisé 
par un coefficient de frottement visqueux a. Etudier son mouvement à l’aide des équations de 
Lagrange. 

Exercice 11 : Etablir l’équation différentielle du mouvement, dans un plan vertical, d’une 

masse ponctuelle m reliée à un point O par une tige de longueur i et de masse négligeable. La 
masse est soumise à une force F(t) qui reste perpendiculaire à la tige lors du mouvement. Les 
forces de frottement de viscosité peuvent être ramenées à une force / = —a v appliquée à la 
masse m dont la vitesse instantanée est v . Le coefficient de frottement visqueux a est supposé 
constant. 
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Chapitre 2 

Oscillations libres des systèmes à un 
degré de liberté 


2.1 Oscillations non amorties 

2.1.1 Oscillateur linéaire 

Un système oscillant à un degré de liberté est habituellement repéré à l’aide d’une coordonnée 
généralisée q qui est l’écart par rapport à la position d’équilibre stable. Le mouvement vibratoire 
est dit linéaire s’il est régi par une équation différentielle harmonique de la forme : 

9 + w oQ = 0 

Cette équation est appelée équation différentielle de l’oscillateur harmonique simple. 


2.1.2 Energie cinétique 

Dans le cas d’un système à un degré de liberté, constitué d’une masse m dont la position est 
repérée par la coordonnée généralisée q, l’énergie cinétique s’écrit : 


1 \dr\ 1 \dfdq] 1 \dr 


T = -mv = -m — = - m — — = - m — q 

2 2 [dt\ 2 [dqdt J 2 [dq\ 

L’énergie cinétique d’un système à un degré de liberté est fonction de q et q . Elle peut s’écrire 

sous la forme : 


T= 2 a(q)q 2 

où a(q) est une fonction de la coordonnée généralisée q, définie dans le cas étudié par : 


dr 

a(q) =m — 


En faisant un développement limité de a(q) au second ordre en q, au voisinage de q = 0 , on 
obtient : 

rr ,. . 1 . . da 1 d 2 a 9 9 

T{q,q) = ~ a{0)+— q +••• Q 


2 dq\ q= ^ 2 dq*\ q=0 

En limitant l’approximation au second ordre, on obtient : 


rr 1 -2 

T = ~a 0 q 


où a 0 est une constante égale à a (0) 
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Oscillations libres des systèmes à un degré de liberté 


2.1.3 Energie potentielle 

Les oscillations se font autour de la position d’équilibre stable q = 0 caractérisée par : 


dU 
dq 

Il est toujours possible , lorsque les écarts par rapport à la position d’équilibre sont faibles, de 
faire un développement en série de Taylor de U(q) au voisinage de la position d’équilibre q = 0. 
En négligeant les puissances de q d’ordre supérieur à deux, on obtient : 



du 

m = m + 


9=0 


1 &U_ 
^ 2 dq 2 


q 2 + 


9=0 


q = 0 correspond à un minimum de U (q) pour lequel 

d 2 U 


dU 

dq 


= 0 et 


q = o 


dq 2 


> 0 


9=0 


Si on choisit l’origine de l’énergie potentielle à cette position d’équilibre (17(0) = 0) , l’énergie 
potentielle U (q) peut s’écrire sous une forme quadratique : 


U(q) \b ü q 2 


, d 2 U 
avec : 6 0 = 


q = 0 


2.1.4 Equation différentielle 

L’équation de Lagrange s’écrit : 


dt 


dC 

dq 


dC . 
“ ~dq ~ J 


Ce qui permet d’obtenir l’équation différentielle de l’oscillateur harmonique simple avec la 
valeur de la pulsation propre ujq : 



6o 

ao 


d 2 u 



a o 


Les oscillations d’un système vibratoire s’effectuent autour d’une position d’équilibre stable. 
Pour des oscillations de faible amplitude autour de la position d’équilibre, tous les mouvements 
vibratoires peuvent être assimilés à des vibrations linéaires et l’énergie potentielle peut alors être 
approximée par une forme quadratique de la coordonnée q, tandis que l’énergie cinétique peut 
être approximée par une forme quadratique en q. 


2.1.5 Résolution de l’équation différentielle de l’oscillateur harmonique simple 

L’équation différentielle de l’oscillateur harmonique simple s’écrit : 

q + u%q = 0 

La solution d’une telle équation est une fonction sinusoïdale du temps 

q(t) = A cos (uiq t + (p) 
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2.2 Oscillations libres des systèmes amortis à un degré de liberté 
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où A représente l’amplitude des oscillations, <p est la phase initiale. 

Il est important de remarquer que la pulsation propre coq ne dépend que des éléments qui 
constituent le système physique étudié (masse, ressort, etc...) tandis que l’amplitude A et la 
phase initiale p sont calculées à partir des conditions initiales : 

q(t = 0) = q 0 
q{t = 0) = q 0 

Enfin l’amplitude des oscillations d’un oscillateur harmonique libre ne dépend pas du temps. 
De telles oscillations sont dites non amorties. 

Il faut néanmoins remarquer qu’au delà d’une certaine amplitude la vibration devient non 
linéaire. Il s’ensuit d’abord une modification de la période des oscillations et ensuite un change- 
ment de la nature du mouvement. 

2.2 Oscillations libres des systèmes amortis à un degré de liberté 

Dans le paragraphe précédent, nous n’avons pas tenu compte de certaines réalités physiques. 
En effet, nous n’avons pas pris en compte les forces de frottement qui sont à l’origine de la perte 
d’énergie mécanique du système sous forme de chaleur. Dans ce paragraphe, nous allons tenir 
compte de ces réalités en nous limitant toutefois au cas simple où les pertes sont dues à des 
frottements visqueux pour lesquels les forces de frottement, qui s’opposent au mouvement, sont 
proportionnelles à la vitesse. 

2.2.1 Equation de Lagrange pour les systèmes dissipatifs 

Rappelons l’équation de Lagrange associée à un système à un degré de liberté dont l’évolution 
au cours du temps se ramène à l’étude de la coordonnée généralisée q 

d_ rac 

dt dq 

Fq représente la composante suivant q de la résultante des forces généralisées qui ne dérivent 
pas d’un potentiel. 

Nous nous intéressons au cas particulier des forces de frottement définies par la force géné- 
ralisée 



F q = f q = -üq 

où (3 est une constante réelle positive. 

L’équation de Lagrange s’écrit alors dans ce cas : 

d_ rac 

dt dq 



2.2.2 Cas particulier des oscillations de faible amplitude 

Nous avons montré dans le chapitre précédent que dans le cas des oscillations de faible 
amplitude, la fonction de Lagrange s’écrivait sous la forme : 


L’équation différentielle du mouvement s’écrit alors : 

a q + bq = —/3g 
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C’est une équation différentielle du second ordre à coefficients constants qui peut se mettre 
sous la forme : 

q -\~ 2 ô q -\- uj q q = 0 

où 6 est un coefficient positif, appelé facteur (ou coefficient) d’amortissement et défini par : 


ojq est la pulsation propre définie par 


5 = 


2 ao 



2.2.3 Résolution de l’équation différentielle 

La solution de l’équation différentielle dépend de la valeur de ô par rapport à uiq : 

- Si 5 > ojq , on dit que le système est suramorti ou apériodique. 

- Si ô = ojq , on dit que l’on a un amortissement critique. 

- Si 5 < üjo , on dit que le système est sous-amorti ou pseudopériodique. 

Cas où le système est suramorti (ô > cuo) 

La solution de l’équation différentielle s’écrit dans ce cas : 

q{t) = A 1 ei~ 5 ~ + A 2 ei~ S+ ^ Z ^\ t 

A\ et A'2 sont des constantes d’intégration définies par les conditions initiales. La figure ci- 
dessous représente q en fonction du temps dans le cas particulier où q( 0 ) = qo et q( 0 ) = 0 . 
q(t) est une fonction qui tend exponentiellement (sans oscillation) vers zéro. 



Régime fortement amorti : variation de q en fonction du temps 

Cas de l’amortissement critique (<5 = cuo) 

La solution générale de l’équation différentielle est de la forme : 

q(t) = (A 1 + A 2 t) e~ 5t 

Dans le cas particulier où q( 0) = qo et q( 0) = 0, 

q{t) = q 0 (l-\-5t) e~ st 

q(t) est encore une fonction qui tend vers zéro sans oscillation lorsque le temps augmente. 
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Amortissement critique : variation de q en fonction du temps 


Cas où le système est sous-amorti ( ô < wo) 

La solution générale de l’équation différentielle est de la forme : 

q(t ) = A e~ 8t cos (c oa t + (/>) 

avec u a = \/ oJq — A 2 ; A et 0 sont deux constantes d’intégration déterminées à partir des 
conditions initiales. Dans le cas particulier où q( 0) = qo et ç(0) = 0, on obtient : 


A = — q 0 

LOA 


= — arctan ( — 



Système faiblement amorti : variation de q en fonction du temps 


2.3 Exercices 

Exercice 1 : Calculer la fréquence des oscillations pour chacun des systèmes suivants dans 

lesquels la masse m est astreinte à un mouvement vertical uniquement : 
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Exercice 2 : Une masse ponctuelle m glisse sans frottement sur une table horizontale. Elle 
est fixée à deux bâtis fixes par deux cordes de masse négligeable tendues horizontalement. En 
supposant que la tension T des cordes reste constante lors du mouvement, calculer la période 
des oscillations pour de faibles amplitudes du mouvement dans la direction x. 



Exercice 3 : Un iceberg de masse volumique pc, assimilable à un parallélépipède régulier et 
homogène de masse M flotte sur de l’eau de masse volumique constante pe ■ Sa surface de base 
est S et sa hauteur est L. 


PG 



On rappelle que la poussée d’Archimède qui s’exerce sur un objet immergé est : Pa = —PeV(J 
où V est le volume immergé et g l’accélération de la pesanteur. 

1. Calculer, à l’équilibre, le volume immergé de l’iceberg en fonction de son volume total. La 
masse volumique de la glace est pc = 900 kg/m 3 ; celle de l’eau est pe = 1000 kg/m 3 . 

2. L’iceberg est écarté d’une distance verticale h par rapport à sa position d’équilibre. Calculer 
la période de ses oscillations quand les frottements sont considérés comme négligeables. 
Faire l’application numérique pour L = 150 m, h = 2 m, g = 9.8 m/s 2 . 


Exercice 4 : Une tige d’acier de constante de torsion C est soudée par son extrémité au centre 
d’un disque homogène de masse M et de rayon R. L’autre extrémité est encastrée dans un bâti 
fixe. Une masse m est soudée au point le plus bas du disque. 



On tourne le disque d’un angle </> o et on le lâche sans vitesse initiale. Déterminer l’expression en 
fonction du temps de l’angle (f>{t) d’écart du système par rapport à sa position d’équilibre. On 
néglige la flexion de la tige d’acier. 

Exercice 5 : Un métronome est schématisé sur la figure ci-dessous. La masse M est soudée à 
l’ extrémité de la tige. La position de la masse m sur la tige peut être réglée. La tige est supposée 
de masse négligeable ; elle est mobile sans frottements autour de O. La masse M étant en bas, 
on l’écarte d’un angle 9q petit et on l’abandonne sans vitesse initiale. 
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1. Quelle(s) condition(s) doit satisfaire le système pour qu’il puisse osciller? 

2. Déterminer l’expression de la période pour des oscillations de faibles amplitudes. 

3. Sachant que M = 80 g, m = 20 g et L = 4 cm, déterminer la distance l pour que la période 
du métronome soit égale à 2 s. 

4. On veut augmenter la période d’oscillation du métronome. Faut-il rapprocher ou éloigner 
la masse m du point O ? 

Exercice 6 : Dans les figures ci-dessous, une tige homogène de masse M et de longueur L 
oscille sans frottement, dans un plan vertical, autour d’un axe fixe perpendiculaire au plan du 
mouvement en O. 

1. Quelle est la déformation du ressort à l’équilibre, sachant qu’à cette position 9 = 0? 

2. Etablir l’équation différentielle du mouvement dans le cas des mouvements de faible am- 
plitude. 

3. A quelle condition le système de la figure (b) peut-il osciller? Quelle est la nature du 
mouvement lorsque cette condition n’est pas satisfaite ? 



Exercice 7 : Quand l’électron d’un atome d’hydrogène, se déplace d’une petite distance x à 

partir de la position d’équilibre, il subit une force de rappel donnée par : 

2 

F = -kx, avec k = 

où r = 0.05 nm correspond au rayon de l’atome. Calculer la pulsation propre loq des os- 
cillations de l’électron. On donne e = 1.6 x 10 -19 C ,m e = 9.1 x 10 -31 kg, en = 8.85 x 
îo- 12 ^- 1 ™- 2 C 2 . 

Exercice 8 : Calculer la période des oscillations d’une particule de charge q et de masse 

m astreinte à se déplacer selon une trajectoire rectiligne entre deux charges égales q fixées en 
x = ±a. 


Exercice 9 : Une particule de masse m se déplace dans un champ de force conservatif avec une 
énergie potentielle donnée par : 

1/ ( r ) = ! ? k ~ æ2 ) p0ur \ x \ <a 
' \ 0 pour \x\ > a 

où a et k sont des constantes. Sachant que a > 0, étudier les types de mouvement possibles selon 
le signe de k. 

Exercice 10 : L’énergie potentielle d’une particule de masse m est 


V(x) 


ex 

x 2 + a 2 


où c et a sont des constantes positives. Représenter graphiquement V en fonction de x. Etudier 
le mouvement des oscillations de faible amplitude au voisinage de la position d’équilibre stable. 
Sachant que cette particule démarre de sa position d’équilibre stable avec une vitesse v, trouver 
les valeurs de v pour lesquelles : 
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1. elle oscille au voisinage de la position d’équilibre ; 

2. elle s’échappe vers +oo ; 

3. elle s’échappe vers — oo. 


Exercice 11 : 

force F(x) : 


Une particule de masse m se déplace dans la région x > 0 sous l’action d’une 


F ( x ) = —moj 2 


x — 


où ui et a sont des constantes. Représenter graphiquement l’énergie potentielle en fonction de 
x. Calculer la période des oscillations de faible amplitude au voisinage de la position d’équilibre 
stable. La particule démarre de cette position avec une vitesse v. Trouver les valeurs de x limitant 
la région des oscillations. Montrer que la période des oscillations est indépendante de v. (Astuce 
pour le calcul de l’intégrale : faire le changement de variable y = x 2 ) 


Exercice 12 : Un bloc de masse 25 kg est monté sur un support en caoutchouc, de masse 
négligeable, qui se comprime de 6.1 cm sous ce poids. Quand le bloc vibre librement, on enregistre 
les positions de la masse après l’avoir déplacé de 5 cm à partir de sa position d’équilibre (voir figure 
ci-dessous). Sachant que le tapis de caoutchouc peut être symbolisé par un ressort de raideur K 
associé à un amortisseur de coefficient de frottement visqueux a , calculer ces coefficients I\ et 
a. 



0,0 0,5 1,0 1,5 2,0 

t(s) 


Exercice 13 : Le système de la figure ci-dessous est constitué d’un cylindre homogène de masse 
M et de rayon R en rotation autour de son axe de révolution fixe (A). Un ül inextensible, de 
masse négligeable, entraîne le cylindre sans glissement sur sa périphérie ; ses deux extrémités 
sont reliées à un bâti fixe (B) par un ressort de raideur K et un amortisseur de coefficient de 
frottement visqueux a. Quelle la valeur critique du coefficient ac ? 



Exercice 14 : Le système mécanique de la figure ci-dessous est constitué d’une tige rectiligne 
AD, homogène, de masse M = 3 kg et de longueur L = 2m. Cette tige peut tourner, dans le 
plan vertical, sans frottement, autour d’un axe horizontal (A) fixe. Les extrémités A et D de 
la tige sont reliées au bâti fixe B 2 par deux amortisseurs identiques de coefficient de frottement 
visqueux a. Le point C, milieu de la tige, est relié au bâti B 1 par un ressort de raideur k. A 
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l’équilibre, la tige est horizontale. Lorsque la tige est écartée de sa position d’équilibre d’un angle 
6*o puis lâchée sans vitesse initiale, elle prend un mouvement oscillatoire amorti de pseudo-période 
1 s. On constate qu’au bout de 5 pseudo-périodes, l’amplitude est égale à 20 % de l’amplitude 
initiale. En déduire la valeur numérique de a puis celle de k. 
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Chapitre 3 

Oscillations forcées des systèmes à un 
degré de liberté 


3.1 Equation différentielle 

Rappelons la forme générale de l’équation de Lagrange pour les systèmes à un degré de 
liberté : 

d raci de dv_ 

dt Ôq dq dq qext 

où F qext est la force généralisée associée à F ext et où la fonction dissipation est V = dq 2 . 

Pour les oscillations de faible amplitude, la fonction de Lagrange pouvait se mettre sous une 
forme quadratique de q et q 

r 1 -2 1 . 2 

F = - 2»o q 

D’où l’équation différentielle du mouvement 

a 0 q + /3q + b 0 q = F qext 

Cette équation peut se mettre sous la forme d’une équation différentielle du second ordre à 
coefficients constants, avec second membre 

q + 2ôq + ujQq = A(t) 


avec 


ù 

2ao 



A(t) 


Fqext 

«o 


3.2 Système masse-ressort-amortisseur 

Considérons l’exemple mécanique de la figure ci-dessous soumis à une force extérieure F(t) 
appliquée à la masse m. 
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Oscillations forcées des systèmes à un degré de liberté 



Système masse-ressort-amortisseur 


Calculons la force généralisée F x conjuguée de la coordonnée x. Pour cela nous pouvons 
utiliser l’une des deux méthodes suivantes : 

- Soit calculer le travail SW de la force F(t) pour une variation Sr de son point d’application 

SW = F ■ ôr = F dx 


On en déduit la x-composante de la force extérieure 

F * = ^ = 

Soit utiliser la définition de la force généralisée 

-* dr 

F - = F ~àï = nt) 


L’équation différentielle du mouvement s’écrit alors 


avec : 


x + 2ôx + üüqX = A(t) 


ô 


a 

2 m 


cuo 



m 


Fit) 

m 


3.3 Solution de l’équation différentielle 

La solution de cette équation différentielle du second ordre est égale à la somme de la solution 
de l’équation sans second membre (ou solution homogène) xn(t) et d’une solution particulière 
de l’équation avec second membre xp(t) : 

x(t) = xn(t) + xp(t) 

Nous avons déjà étudié l’équation sans second membre xn(t) et nous savons que cette solution 
contient dans tous les cas le terme exponentiel e~ St . Après un intervalle de temps t supérieur 
à 3/èou 4/é, le terme e~ St devient très petit et la solution homogène est alors pratiquement 
nulle. Il ne subsistera que la solution particulière de l’équation avec second membre. L’intervalle 
de temps pendant lequel la solution homogène est non négligeable est appelé le régime transitoire. 
A la fin de ce régime transitoire commence l’intervalle de temps pour lequel la solution homogène 
est quasi-nulle et pour lequel la solution x(t) — x p (t) ; ce régime est appelé régime permanent ou 
stationnaire. 
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3.3.1 Cas particulier où A(t) = A 0 cos(fit) 

Calcul de la solution permanente à l’aide de la méthode des nombres complexes 

Pour t suffisamment grand, nous pouvons considérer que la solution transitoire s’est annulée 
et que la solution x(t) s’identifie alors avec la solution particulière : x(t) — xp(t). Par commodité 
de notation l’indice p est sous-entendu dans ce qui suit. La méthode des nombres complexes 
permet de calculer aisément la solution stationnaire. 

Soit le déplacement complexe représenté par le nombre complexe X = Je“, avec X_ = 
Xç>e l,f . Nous pouvons considérer, en outre, que A(t.) = Aq cos (fit) constitue la partie réelle du 
nombre complexe A = L’équation différentielle se transforme en une simple équation 

algébrique en fonction de l’amplitude complexe X : 


[(wg - fi 2 ) +i25fl\X = A 0 


dont la solution est : 


X = 



Ap 

fi 2 ) +i2ô fi 


D’où l’on tire l’amplitude Xç> et la phase cp : 


X 0 = 


Aq 

(uj 2 - fi 2 ) 2 + 4Ô 2 fi 2 


2 S fi 

p = — arctan — - 

cjg — ii 2 


Etude des variations de l’amplitude et de la phase en fonction de la pulsation de 
l’excitation 


Le maximum de l’amplitude est obtenu pour la valeur de fi qui annule 


dX p 
dfi 


Il existe un maximum à la pulsation fip = y/l ug — 2 ô 2 seulement si l’amortissement est 
suffisamment faible pour que ô < ujq//2. A cette pulsation, appelée pulsation de résonance, on 
dit que le sj^stème entre en résonance et l’amplitude Xq est maximale ; elle vaut : 


V - A ° 

° maX 2 ôy/l^P 

La figure représentant les variations de Xo en fonction de la pulsation d’excitation fi est 
appelée courbe de résonance en amplitude. On remarque qu’à la pulsation uiq, le déphasage p 

TT 

est égal à — — , et qu’à la résonance ip = — arctan 
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Amplitude X 0 en fonction de fl. Déphasage p en fonction de fl 


Etude de la résonance pour les faibles amortissements 

Dans le cas des faibles amortissements ( ô « wo), la fréquence de résonance est très peu 
différente de la pulsation propre, Qr — ujp . Dans ce cas, l’amplitude de vibration à la résonance 
A'omax est égale à : 


y _ ^0 

^ 0 max — 0 e 

ZÔCüo 

Pour les faibles amortissements, Xo max est donc inversement proportionnel à ô. 


Etude de la vitesse 


En notation complexe, la vitesse s’écrit : 

V(t) = — = iÜX = Xe iÇit 
dt 

où l’amplitude complexe de la vitesse est définie par 

X = ÎÜX = 

— — (u%-n 2 )+i 2ÔQ 

L’étude des variations de l’amplitude de la vitesse en fonction de la pulsation d’excitation 
montre que, quelle que soit la valeur de ô, la résonance en vitesse est obtenue pour fl = luq (voir 
figure ci-dessous). La valeur maximale de l’amplitude de la vitesse vaut dans ce cas : 


X max — X(cco) 


Ap 

25 



Courbe de résonance de la vitesse 



Déphasage ip de la vitesse en fonction de 
fl. 
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Bilan énergétique 

Soit Pp(t) la puissance instantanée fournie par la force extérieure F(t) au système. En régime 
permanent, on obtient : 

Pp{t) = F(t ) x(t) = FqXq cos (fît) cos (fît + V 1 ) 

Soit < Pp > la valeur moyenne sur une période de Pp(t) : 

< Pf >= ^F 0 X 0 cos (■!/’) 

En tenant compte de l’expression de Xq en fonction de Fq, on obtient : 

1 ■ 9 

< Pf >= 2 aX o 

Comparons cette valeur à la valeur moyenne < Pd > de la puissance dissipée par les forces 
de frottement de viscosité. La valeur instantanée de cette puissance dissipée s’écrit : 

Pn(t) = ax 2 = aXfj cos 2 (fît + tp) 

D’où l’on tire la valeur moyenne sur une période : 

1 ■ n 

< Pd >= ÿ a Xo 

L’étude des variations de la valeur moyenne de la puissance < P >=< Pp >=< Pd > en 
fonction de la pulsation d’excitation montre que la valeur maximale de la puissance moyenne est 
obtenue pour fî = uq quelle que soit la valeur de S. La valeur maximale de la puissance moyenne 
dissipée ou fournie vaut dans ce cas 


< P > 


max — 


2a 


La figure ci-dessous représente les variations, en fonction de fî, de la puissance moyenne 
dissipée par les forces de frottements (ou de manière équivalent la puissance moyenne fournie 
par la force extérieure ) . 



Courbe de résonance pour la puissance 
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Bande passante 

On définit par bande passante, la bande des pulsations autour de fî = ojq pour lesquelles 
< P >>< P >max /2. Les deux pulsations Qi et fl 2 , situées de part et d’autre de la pulsation 
cüo et pour lesquelles < P >=< P > max /2, sont appelées pulsations de coupure. La bande 
passante B s’écrit : 


b = n 2 -n i 


Le calcul de B consiste à rechercher les deux pulsations pour lesquelles < P >=< P > max /2. 
On obtient l’expression de la bande passante B : 

B = ü 2 — fii = 26 


Coefficient de qualité d’un oscillateur 

Le coefficient de qualité d’un oscillateur est défini par le rapport de la pulsation propre ojq à 
la largeur de bande B : 


3.3.2 Cas d’une excitation périodique 

Nous avons étudié dans le paragraphe précédent la réponse d’un système vibratoire à une 
excitation sinusoïdale dite excitation harmonique. En pratique, les excitations mécaniques ne 
sont pas toujours parfaitement sinusoïdales ; elles sont souvent périodiques. En considérant le 
cas d’excitations périodiques, nous procéderons à une généralisation du cas harmonique. 

Soit une excitation périodique appliquée à un système amorti à un degré de liberté. L’équation 
différentielle qui régit ce système s’écrit : 

q + 25q + ulq = A(t ) 

La fonction A(t) étant périodique, de période T, son développement de Fourier s’écrit : 


OO 

A(t) = — + ^2 a n cos(nut) + b n sin(ntcf) 

n= 1 

L’équation différentielle s’écrit alors : 


q + 28 q + ojQq 


a o 

T 


OO 

+ ^2 a n cos(nut) + b n sin (nut) 

n= 1 


La réponse permanente (ou stationnaire ) qui s’identifie avec la solution particulière, pour t 
suffisamment élevé, peut alors être calculée pour chacune des composantes de l’excitation : üq/2, 
a n cos(nut) et b n sm(nujt). On obtient alors par superposition : 


OO 

u 1 


a n cos (a ] n t + ipn) + bn sin (a > n t + ip r , 
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3.4 Impédance mécanique 

3.4.1 Définition 

Considérons un système mécanique soumis à une force sinusoïdale F (t) = Focos(Pt). En 
régime permanent, le point d’application de cette force se déplace avec une vitesse v ( t ) = 
Vo cos (Pt + (f>) . On appelle impédance mécanique d’entrée du système mécanique, le rapport 
des amplitudes complexes de la force F et de la vitesse v 


Z.E ~ 


F 

V 


3.4.2 Impédances mécaniques 
Amortisseur 

Dans le cas d’un amortisseur, la force appliquée est reliée à la vitesse par 

F = av 

On en déduit l’impédance complexe d’un amortisseur 

Z_ a = a 


Masse 

Dans le cas d’une masse, la relation fondamentale de la dynamique s’écrit 

. dv 

F = m— 
dt 


On en déduit l’impédance complexe d’une masse 

Z_ m = irriil = mQ e* 2 


Ressort 

Dans le cas d’un ressort de raideur k, la force appliquée / appliquée au ressort s’exprime en 
fonction de l’allongement par 

f = kx 

On en déduit l’impédance complexe d’un ressort 




k 

n 


e -*2 


3.4.3 Puissance 

La valeur moyenne, sur une période, de la puissance fournie est 

< Pf >= -^FqXq cos (<f>) = - Re ( Zg) Xf } 
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3.4.4 Applications 
Système mécanique résonant 

Soit un système mécanique constitué d’un ressort de raideur k, d’un amortisseur de coefficient 
de frottement visqueux a et d’une masse m soumise à une force sinusoïdale F (t) = Fq cos (f U). 
L’impédance d’entrée de ce système est 


Z f — oz i 


mil - ^ 



le module de l’impédance est Ze = a. Lorsque la pulsation 


A la résonance 
Q — > oo, l’impédance Z_e — imil. 



Module de l’impédance d’entrée. 



Système antirésonant 

Considérons un système mécanique constitué d’un ressort de raideur k dont une extrémité est 
reliée à une masse m et dont l’autre est soumise à une force sinusoïdale F (t). Soit x le déplacement 
de la masse m et soit y le déplacement du point d’application de la force F(t). Pour calculer 
l’impédance d’entrée de ce sj^stème, nous devons d’abord écrire les équations différentielles du 
mouvement : 

mx = k (x — y ) 

F = k (x — y) 

En utilisant la notation complexe, on obtient l’impédance d’entrée : 



km 




mü - ^ 


La pulsation d’antirésonance est u o = \ — . Lorsque ü, = ujq, 
que le module de l’impédance est oo. Lorsque la pulsation fî — >• oo, 


la vitesse Ÿ est 
l’impédance Ze 


nulle tandis 
-> 0 . 
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3.5 Exercices 

Exercice 1 : Un disque circulaire homogène, de masse M, de rayon R, peut osciller sans 

frottements autour de son axe horizontal O. Deux masses m\ et rri 2 sont soudées aux extrémités 
d’une tige de masse négligeable liée rigidement au disque et passant par O. Les distances de m\ 
et ??t -2 au centre sont notées respectivement l\ et £ 2 • Un ressort vertical, de constante de raideur 
K a une extrémité fixe et l’autre est reliée au disque en un point A situé à une distance a de O. 
En position d’équilibre la tige est verticale avec m\ en bas et le point A est au même niveau que 
le centre O. Le disque subit un frottement visqueux de coefficient a au point B. La masse m\ est 
soumise à une force F(t) = Fq cos(flf) perpendiculaire à la tige. Valeurs numériques : M = 1 kg, 
777-1 = m 2 = 0.1 kg, K = 16 N/m, R = 20 cm, £\ = 50cm, £2 = 25cm, a = 10 cm, g = 10 m/s 2 , 
a = 7.25 x 10 -2 kg/s. 



1. Etablir l’équation différentielle du mouvement. 

2. Trouver sa solution en régime permanent. 

3. Calculer le facteur de qualité Q du système. 

4. Déterminer la valeur de Fq pour qu’à la résonance l’amplitude maximale soit égale à 7r/30 
rad. 


Exercice 2 : (suite de l’exercice n°14 du chapitre précédent) Le bâti B\ est maintenant animé 
d’un mouvement vertical sinusoïdal donné par : s(t ) = Sqcos(ÇH') où Sq = 1 cm. 
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1. Montrer que l’équation différentielle qui régit le mouvement du système peut s’écrire : 

6 + 2 ô 9 + ujq 9 = Aq cos(f U) 

On précisera de manière explicite le terme Aq. Calculer sa valeur numérique. 

2. Quelle est l’expression de la solution 6(t) lorsque le régime permanent est établi? Véri- 
fier que le système est très faiblement amorti ; en déduire la fréquence de résonance et 
l’amplitude de 6(t) à la résonance. 

3. Quelle est, à la résonance, l’amplitude de la force Ft transmise au sol par chaque amortis- 
seur ? 


Exercice 3 : Le dispositif mécanique ci-dessous représente le schéma de principe d’un appareil 
de mesure de vibrations. La masse m est liée par deux ressorts et un amortisseur de coefficient de 
frottement visqueux a à un support rigidement lié au système mécanique dont on veut étudier les 
vibrations. Le mouvement du support est repéré par s(t) tandis que le mouvement de la masse 
est repéré par x(t). On étudie des vibrations sinusoïdales de la forme s(t) = So cos(fii). L’origine 
est prise à la position d’équilibre. 



1. Etablir l’équation du mouvement de la masse m en fonction de la coordonnée relative 
y{t) = x(t) - s(t). 

2. Déterminer la solution stationnaire y(t). 

3. Dans le cas de ressorts de faible raideur, la pulsation propre uq est petite devant la pul- 
sation D. Donner dans ce cas l’expression de y(t). Montrer que l’on peut ainsi déterminer 
facilement l’amplitude So de la vibration (on a réalisé ainsi un vibromètre). 

4. Lorsque la raideur des ressorts est élevée, la pulsation propre luq est grande devant la pulsa- 
tion Q des vibrations. Montrer , que dans ce cas on peut déterminer facilement l’accélération 
du support (on a ainsi réalisé un accéléromètre). 

Exercice 4 : Un véhicule roulant est un système complexe à plusieurs degrés de liberté. La 
figure ci-dessous peut être considérée comme une première approximation d’un véhicule qui se 
déplace sur une route ondulée décrite par le profil yi(t). 



Dans ce modèle simplifié, on suppose que : 

- La raideur élastique des pneus est infinie, c’est-à-dire que les ondulations de la route sont 
intégralement transmises à la suspension du véhicule. 
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- Les roues ne décollent pas de la chaussée. 

- On s’intéresse uniquement au déplacement vertical y(t) du véhicule dans le plan de la 
figure. 

- On se place dans le cas simple où le véhicule se déplace horizontalement à une vitesse 
constante v sur une route à profil sinusoïdal yi(x) = Y\ sin(27rx/A). 

1. Etablir l’équation différentielle qui régit les variations au cours du temps de la coordonnée 
y du véhicule. 

2. En déduire l’amplitude Y du mouvement du véhicule dans le sens vertical . 

3. Application numérique m = 350 kg, k = 350 kN/m, v = 100 km/h, A = 5m, Y\ = 20 cm ; 

(a) pour a = 2000 N s/m, 

(b) pour a = 200 N s/m. 


Exercice 5 : Les machines tournantes (moteurs électriques, turbines, machines à laver, etc...) 
peuvent être le siège de vibrations importantes car très souvent le centre de masse ne coïncide 
pas avec l’axe de rotation. Pour limiter ces vibrations on utilise des supports antivibratoires 
constitués généralement de caoutchouc renforcé. En raison de leurs propriétés mécaniques ces 
supports peuvent être modélisés par un amortisseur en parallèle avec un ressort. 

On se propose d’étudier à titre d’exemple le cas d’une machine à laver le linge (figure 1 ci- 
dessous). Soit M la masse de cette machine. La partie tournante est constituée d’un tambour de 
rayon e tournant à une vitesse angulaire constante f 2. On considère que la masse tournante est 
constituée par le linge de masse m. Pour des raisons de simplicité, on suppose que le lave-linge 
ne peut effectuer que des mouvement verticaux repérés par la coordonnée y. 


M M+m 



Figure 1 


Figure 2 


1. Etablir l’équation différentielle du mouvement pour la coordonnée y. 

2. Montrer qu’un tel dispositif est équivalent au schéma simplifié de la figure 2 ci-dessus ; 
donner l’expression de F eq . 

3. Dans l’hypothèse des faibles amortissements ( 5 « iuo ), tracer et commenter le graphe de 
l’amplitude Y du déplacement vertical du lave-linge en fonction de la vitesse de rotation. 

4. Calculer l’amplitude de la force transmise au sol à la résonance. 


H. Djelouah 





28 


Oscillations forcées des systèmes à un degré de liberté 


H. Djelouah 



Chapitre 4 

Oscillations libres des systèmes à deux 
degrés de liberté 

4.1 Introduction 

Les systèmes qui nécessitent deux coordonnées indépendantes pour spécifier leurs positions 
sont appelés systèmes à deux degrés de liberté. 

Exemples 



- Figure 1 Si les masses m\ et m 2 sont astreintes à se déplacer verticalement, 2 coordonnées 
x\ et X 2 sont nécessaires pour spécifier la position de chaque masse à chaque instant. 

- Figure 2 Si la masse M est astreinte à se déplacer dans un plan vertical, deux coordonnées 
sont nécessaires pour spécifier la configuration du système. L’une de ces coordonnées peut 
être le déplacement x qui correspond à la translation verticale de la masse. L’autre coor- 
donnée peut être le déplacement angulaire 9 pour tenir compte de la rotation de la masse. 
Ces deux coordonnées sont indépendantes l’une de l’autre. 

- Figure 3 Dans le cas du double pendule, deux coordonnées sont nécessaires pour spécifier 
la position des masses m\ et 777 , 2 • Plusieurs choix sont pourtant possibles, en effet on peut 
choisir (xi,x 2 ) ou (yi,y 2 ) ou (9i,9- 2 ). 

Il est possible de spécifier la configuration d’un système à l’aide de plusieurs ensembles de 
coordonnées indépendantes ; un ensemble quelconque de ces coordonnées est appelé coordonnées 
généralisées. Il y a autant d’équations de Lagrange que de degrés de liberté ou de coordonnées 
généralisées. Pour l’étude des systèmes à deux degrés de liberté, il est nécessaire d’écrire deux 
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équations différentielles du mouvement que l’on peut obtenir à partir des équations de Lagrange 


' d 

de 

dt 

Ôqi 

d 

'dC 

dt 

dq 2 _ 


()£ 

dqi 


= 0 


dC 

dq2 


4.2 Systèmes à deux degrés de liberté 

4.2.1 Système masses- ressorts en translation 

m i 7712 



Système masses-ressorts en translation 

Considérons le système ci-dessus, constitué de deux masses m i et m 2 reliées respectivement 
par deux ressorts de raideur k\ et k 2 à deux bâtis fixes. Les deux masses sont reliées par un 
ressort de raideur K. Ce ressort est appelé ressort de couplage. 


Equations différentielles du mouvement 

Les équations du mouvement pour ce système à deux degrés de liberté peuvent être obte- 
nues à partir des équations de Lagrange pour chaque coordonnée x\ (t) et X2 (t). Soit T et U 
respectivement l’énergie cinétique et l’énergie potentielle : 

T = \ m\ x'I + \ iri2 x\ 

U = \k\ x\ + \K (aq - X2) 2 + \k 2 x\ 


U = \ (k\ + K) x\ + \ ( k 2 + I< ) x\ - Kx\X2 
Le lagrangien C = T — U s’écrit alors 

C = ^ mi x\ + ^ m 2 x\ - ^ (h + K ) x\ - ^ (k 2 + I<) x\ + I\xix 2 
Les éc|uations de Lagrange s’écrivent 

' <1 
dt 

< 

A. 

dt 

D’où le système d’équations différentielles du mouvement 


dC 

dxi 

dC 

dx 2 


dC 

dx\ 

dC 

dx 2 


= 0 


= 0 


m\X\ + (ki + I\) xi — Kx 2 = 0 
VI2X2 + (k2 + K) X2 - Kx 1 = 0 


Les termes —Kx 2 et —Kx 1 qui apparaissent respectivement dans la première et la seconde 
équation sont appelés termes de couplage, et les deux équations différentielles sont dites couplées. 
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Résolution des équations différentielles 

Recherchons une solution particulière de la forme : 

aq (t) = A\ cos (uit + 4>) 

X2 (t) = A 2 COS (ait + (j)) 

où Ai, A2 et (j) sont des constantes et ex l’une des pulsations propres du système. La substi- 
tution de aq et X2 dans le système d’équations différentielles donne 

( [ki + K — mi ex 2 ] Ai — K A 2 = 0 

[ - K Al + [k 2 + K - 7772 ex 2 ] A 2 = 0 

Ce qui constitue un système d’équations linéaires homogènes dont les inconnues sont Ai et 
A2 . Ce système admet une solution non identiquement nulle seulement si le déterminant A (ex) 
des coefficients de Ai et A2 est égal à zéro. 

A / \ [h + I< - 7771 W 2 ] - K 

M -K [k 2 + K - m 2 Lü 2 ] 

Le déterminant A(ex) est appelé déterminant caractéristique. L’équation A(ex) = 0 est appelée 
l’équation caractéristique ou équation aux pulsations propres. Elle s’écrit 

[L’I + K — 7771 CX 2 ] \k2 + K — 7772 tx 2 ] — K 2 = 0 


ou encore 


4 2 

ex — ex 


ki + K k 2 + K 


? 77 l 


777 2 


kik 2 + kiK + k 2 K _ o 

777 1 7772 


Cette équation est une équation quadratique en ex qui admet deux solutions réelles positives 
exi et ex2 appelées les pulsations propres du système. 

Cet exemple montre qu’il y a en général deux pulsations propres dans un système à deux 
degrés de liberté. Chacune des coordonnées, aq et X2, possède deux composantes harmoniques 
de pulsations exi et ui 2 


xi = An cos(cxi t + 4 >i) + A12 cos(cx 2 i + (f> 2 ) 
x 2 = A 2 1 cos(cxi t + 4 >i) + A 2 2 cos(t 0 2 t + (f> 2 ') 

où An, A12, A21, A22, 4 >i et 02 sont des constantes. Le terme de plus basse fréquence corres- 
pondant à la pulsation exi est appelé le fondamental. L’autre terme, de pulsation 7x2, est appelé 
harmonique. 

Les doubles indices sont utilisés pour les amplitudes des différentes composantes harmoniques ; 
le premier indice se réfère à la coordonnée et le second à la pulsation. Par exemple An est 
l’amplitude de aq(t) à la pulsation 1x2- 

Lorsque An = A22 = 0 , aq et aq correspondant à la première solution particulière sont des 
fonctions sinusoïdales, en phase, de pulsation exi ; on dit que le système oscille dans le premier 
mode. Dans ce cas 


aq = An cos(cxit + </q) 


x 2 = A 2 i cos(oqf + </q) 
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Lorsque An = A21 = 0, aqet X2 correspondant à la seconde solution particulière et sont des 
fonctions sinusoïdales, en opposition de phase, de pulsation lu 2 ; on dit que le système oscille dans 
le second mode. Dans ce cas 


xi = An cos(uj 2 t + ) 


X2 = A22 COs(c U 2 t + $2) 


Etudions les particularités de ces deux solutions particulières : 
- La première solution particulière s’écrit : 

x\ = An cos(cJit + 4> 1 ) 


X2 = A 2 i cos(cji t + 4 >i) 

x\ et X2 doivent vérifier le système d’équations différentielles, ce qui donne 

( [ki + K - miujj] An - K A 2 \ = 0 

\ - I< An + \k2 + K - m 2 ujf\ A 2 i = 0 

Ces deux équations permettent d’obtenir le rapport des amplitudes dans le premier mode 
ou fondamental 

A 21 ki + K — m 1 ujf K 

^ 1 .A 11 K k 2 + K - m 2 ijj\ 

- La seconde solution particulière s’écrit : 

Xl = A\2 C0S(CJ2 1 + 02) 

X2 = A 22 C0S(CJ2 1 + 0 2 ) 

x\ et X2 doivent vérifier le système d’équations différentielles, ce qui donne 
( [ki + K - rn\ ix An - K A 22 = 0 

[ - I< A 12 + [k 2 + K - m 2 w|] M 2 = 0 

Ces deux équations permettent d’obtenir le rapport des amplitudes dans le second mode 
ou harmonique 

A22 ki + K — m 1 lu I K 

^ 2 A12 K k 2 + K - m 2 LU 2 

- La solution générale (xi,x 2 ) est une combinaison linéaire de ces deux solutions particu- 
lières. xi et X2 s’écrivent alors 

xi = An cos (oui t + 4 >i) + A12 cos (cu 2 t + fc) 


X2 = n 1 An cos (cuit + 4 >i) + /i 2 A i2 cos (cu 2 t + ^2) 

où An, A 12, 4 >i et 4> 2 sont des constantes d’intégration dont les valeurs sont fixées par les 
conditions initiales. 
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4.2.2 Cas particulier de deux oscillateurs identiques 
Calcul des constantes d’intégration 

Considérons le cas particulier de deux oscillateurs identiques tels que m\ = m 2 = m et 
ki = k-2 = k. Dans ce cas les pulsations propres sont respectivement égales à 



, , / k + 2 K , , /ï i 2 K 

w 2 - yj - w i y 1 + X 

Les rapports d’amplitudes correspondant à ces pulsations sont respectivement fi\ = +1 et 
M 2 = -b 

Soit xio,X2o,xio et x 2 o les valeurs initiales respectives de x’i,x 2 ,xi et x 2 ■ Tenant compte de 
ces conditions initiales, on obtient le système d’équations suivant qui permet de déterminer les 
constantes d’intégration An, A12, 4 >i et <f>2 


An cos(<j 4 ) + A12 cos(</> 2 ) = aqo 
An cos((^i) - A V2 cos(4>2) = x 20 
-uji An sin(çi)i) - w 2 4 i2 sin( 0 2 ) = x w 
-oj\ An sin( 0 i) + 0J2 An sin( 0 2 ) = x 20 

Les solutions de ce système d’équations sont 


ou encore 


æio + X20 , A æio — x 20 

= ; et A 12 = ; 

2 cos(^i) 2 cos (fo) 


Au = 


æio + X20 
2 sin(^i) 


et 


4-12 = 


.t 2 q — æio 
2 0J2 sin(</> 2 ) 


1 . Considérons le cas particulier suivant xio = x 2 o = xq et ±10 = x 2 o = 0 ; on obtient dans 
ce cas (f>i — <f>2 — 0 , 4 i 2 = Oet An = xq ; d’où 


x\ = xo cos(cjit) 

X2 = xq cos(cjit) 

Pour ces conditions initiales particulières, les deux masses oscillent en phase à la même 
pulsation On dit que le système oscille dans le mode fondamental. 



Oscillations dans le mode fondamental 
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2. Considérons un autre cas particulier pour lequel xio = — X20 = xq et iqo = X20 = 0 . On 
obtient dans ce cas 0i = 02 = 0, An = 0 et A12 = xq ; d’où 

X\ = XQ COs(u> 2 t) 

X 2 = — Xq C0S(W2 t) 


On dit que le système oscille dans le second mode car les deux masses oscillent en opposition 
de phase avec la même pulsation ui2- 



Oscillations dans le mode harmonique 


3. Considérons enfin le cas particulier suivant xio = xo, X20 = 0 et x\q = X20 = 0; d’où 
</>i = 4>2 = 0, An = A 12 = xq/2. Les solutions s’écrivent alors sous la forme 


xi(t) = ^ cos ( w i t) + ^2 cos ( w 2 t) 

X2 (t) = ^ cos (üji t) — ^ cos (u>2 t ) 

Les solutions ne sont plus des fonctions purement sinusoïdales du temps mais des combi- 
naisons linéaires de deux fonctions sinusoïdales de pulsations respectives uj\ et 1x2- x\ et X2 
peuvent s’écrire sous la forme 


X\ ( t ) = xo cos 
X2 ( t ) = xq sin 



cos 


sin 


0J2 + Wl 
2 

0J2 + OJl 


La figure suivante représente le résultat obtenu dans le cas où uj\ est très différent de 012 
(c’est-à -dire si K » k). 



Oscillations dans le cas des conditions initiales : xio = xq, X20 = 0 et xio = X20 = 0 
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Si wi est peu différent de uj-2 (c’est-à -dire si K « k), on observe un phénomène de 
battement (voir figure ci-dessous). 



Coordonnées principales 

Considérons les coordonnées p\ et p2 obtenues à partir des coordonnées x\ et xq par les 
relations 


Pi 


Xï + XQ 
2 


P2 = 

Tenant compte des expressions de x\ et X2 et des valeurs particulières de p\ et p2 pour 
l’exemple étudié, on obtient 


Pi = cos (uit) 

P 2 = y cos (u 2 t) 

On remarque que, quelles que soient les conditions initiales, p\ et P2 sont des fonctions 
purement sinusoïdales du temps de pulsations respectives et W2- Ces coordonnées particulières 
sont appelées coordonnées principales. On peut vérifier que le système d’équations différentielles 
qui régit le mouvement du système considéré s’écrit sous la forme de deux équations découplées 

pi + uj pi = 0 

P 2 + u 2 p 2 = 0 


Les relations inverses suivantes 


xi=pi+ p 2 
X2=Pl~ P-2 


permettent d’obtenir les coordonnées aq et X2 à partir des coordonnées principales pi et p2- 
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4.2.3 Pendules couplés 

Considérons le cas de deux pendules simples identiques couplés par un ressort de raideur K 
et qui effectuent des oscillations de faible amplitude repérées par les angles 9\ et 9 2 . 



Pendules couplés 


Etablissons tout d’abord les équations différentielles du mouvement dans le cas des oscillations 
de faible amplitude. Il est aisé de montrer que l’énergie cinétique et l’énergie potentielle s’écrivent 
sous les formes quadratiques suivantes 

T = \ml 2 9\ + |m/ 2 9 2 

U = \ [Kl 2 + mgl ] 0\ + \ [Kl 2 + mgl ] Q\ - Kl 2 9i9 2 

On remarque la présence du terme de couplage -Kl 2 6 i 6-2 dans l’expression de l’énergie po- 
tentielle. Comme dans l’exemple précédent, on dit que le couplage est élastique. Si le terme de 
couplage n’existe que dans l’expression de l’énergie cinétique, on dit que le couplage est de type 
inertiel. 

Les équations de Lagrange permettent d’obtenir les équations différentielles du mouvement 

ml 2 ë i + [Kl 2 + mgl] 9 X - Kl 2 6 2 = 0 
—Kl 2 9i + ml 2 92 + [Kl 2 + mgl] 92 = 0 

En l’absence d’amortissement une solution particulière de ce système d’équations différen- 
tielles serait 


9\(t) = Ai cos (ut + (/>) 

6*2 (t) = A 2 cos (u>t + 4>) 

Ces deux expressions doivent satisfaire le système d’équations différentielles, d’où 

[Kl 2 + mgl — ml 2 uj 2 ] Ai — I\l 2 A 2 = 0 
—Kl 2 Ai + [Kl 2 + mgl — ml 2 uj 2 ] A 2 = 0 

H. Djelouah 




4.3 Exercices 


37 


Ce système d’équations admet des solutions non nulles seulement si uj est solution de l’équa- 
tion aux fréquences 

[Kl 2 + mgl - ml 2 J 2 ] 2 - I< 2 l A = 0 
D’où l’on tire l’expression des pulsations propres uq et uq 

[g , 

uq = W y et U 2 = 

La solution du système d’équations différentielles est donc 

8 1 = An cos(uq t + </>i) + j4i2 cos (uqi + <fo) 

8-2 = A 2 1 COS (uqt + </>l) + ^22 COS (w 2 f + <t> 2 ) 

Pour calculer les rapports des amplitudes dans les modes, on suppose que le système oscille 
soit dans le premier mode soit dans le second mode. Dans le premier mode, on obtient le système 

[Kl 2 + mgl — m,l 2 üjf\ — Kl 2 m = 0 

—Kl 2 + [Kl 2 + mgl — ml 2 uf\ = 0 
Dans le second mode, on obtient 

[Kl 2 + mgl — ml 2 u 2 ] — Kl 2 g 2 = 0 
—Kl 2 + [Kl 2 + mgl — ml 2 uf\ g 2 = 0 

Tenant compte des expressions de uq et W 2 on obtient les valeurs du rapport des amplitudes 
dans les modes g \ = +1 et g 2 = — 1. Les solutions du système d’équations différentielles s’écrivent 
alors 

81 = An cos (uq t + 4 >i) + Ai 2 cos + ^ 2 ) 

82 = An cos (uq t + </>i) - ^4 i 2 cos (u 2 t + c/> 2 ) 

4.3 Exercices 

Exercice 1 : Soit le système mécanique représenté par la figure ci-contre, composé de deux 
oscillateurs linéaires (m, k) couplés par un ressort de raideur k' . 



1. Ecrire le lagrangien du système. 

2. (a) Mettre ce Lagrangien sous la forme : 

L = -m [(xf + — u>q (xf + x| — 2Cxix 2 y\ 

Donner les expressions de lüq et C (coefficient de couplage). 
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(b) En déduire les équations du mouvement. 

3. (a) Déterminer les pulsations propres du système. 

(b) Sachant que Tq = 2-k/ujo = 0, 5 s et C = 0,3, calculer les valeurs numériques des 
périodes propres. 

4. (a) Le coefficient de couplage C étant faible, donner les solutions x\ (t) et X 2 (t) avec les 

conditions initiales suivantes : à t = 0 s, x\ (0) = Xq , x \ (0) = 0 et X 2 (0) = 0 ,X 2 (0). 
Quel phénomène physique observe-t-on ? 

(b) Tracer l’allure des courbes représentatives de x\(t) et X 2 (t). 

Exercice 2 : Soit le système mécanique représenté figure ci-dessous. Les variables x\ (t) et X 2 (t) 
représentent les déplacements horizontaux (à partir de l’équilibre) des masses m \ et m ,2 dans le 
cas des petites oscillations. La tige de longueur L est de masse négligeable. 



X 1 *2 

On se place dans le cas où : K\ = K 2 = K :i = k et m\ = m 2 = m. On posera : 

2 / 5k g\ 2k 

Wo= \4^ + L) = m 

1. Calculer les pulsations propres 

2. Déterminer les rapports des amplitudes de chacun des modes. 

3. En déduire l’expression de x\ (t) et X 2 (t). 

4. Donner les solutions de aq(t) et X 2 (t) si x\ (0) = xo,xi (0) = 0 et X 2 (0) = 0,X2 (0) = 0 

Exercice 3 : Soit le sj^stème mécanique suivant comprenant entre autres une barre horizontale 
de masse négligeable et qui peut pivoter sans frottement autour d’un axe passant par son milieu. 



On prendra M = 2m et k\ = k 2 = k% = k 

1. Etablir les équations régissant les petites oscillations. 

2. Trouver les pulsations propres et les rapports des amplitudes pour les différents modes. 

3. Ecrire les solutions générales y(t.) et 9(t). 
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Exercice 4 : Dans la figure ci-dessous, M et R représentent respectivement la masse et le 

rayon de la poulie. x\ et X 2 représentent les écarts des deux masses par rapport à leur position 
d’équilibre. 



On prend : M = 2(m2 — m{) avec rri 2 = m, et ko = k\ = = k. 

1. Ecrire le Lagrangien du système. 

2. Déterminer les pulsations propres et le rapport des amplitudes de chacun des modes en 
fonction de m et k. 


Exercice 5 : Sur la figure ci-dessous, nous avons schématisé un véhicule avec sa suspension 
(sans amortisseurs). Nous supposons que les ressorts restent verticaux. La masse du véhicule est 
m et son moment d’inertie par rapport à un axe horizontal D passant par le centre de gravité G 
et perpendiculaire au plan de la figure est Jq- Le déplacement du centre de gravité par rapport 
à l’équilibre est repéré par x (pompage). L’angle 6 (tangage) que fait le châssis avec le sol, par 
rotation autour de D, sera supposé petit. L’inclinaison sur les côtés (roulis) est supposée nulle. 



On donne les valeurs suivantes : 

- masse du véhicule m = 1000 kg, 

- distance entre l’axe avant et G : L\ = 1 m, 

- distance entre l’axe arrière et G : L 2 = 1.5 m, 

- constante de raideur du ressort avant : k\ = 18 kN/m, 

- constante de raideur du ressort arrière : k 2 = 18 kN/m, 

- moment d’inertie du véhicule : Jq = mr 2 ; r = 0.9 m. 

1. Déterminer les pulsations propres du système ainsi que le rapport des amplitudes dans 
chacun des modes. 

2. Ecrire les solutions x(t) et 6(t). 

3. (a) Quelle condition doit être réalisée si l’on désire avoir un découplage entre x et 9 ? 

(b) Quelles sont alors les fréquences propres de pompage fp et de tangage fp ? 
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Chapitre 5 

Oscillations forcées des systèmes à deux 
degrés de liberté 


5.1 Equations de Lagrange 

Soit un système à deux degrés de liberté, soumis à des forces qui dérivent d’un potentiel, à 
des forces de frottement de viscosité et à des forces extérieures. Si les coordonnées généralisées 
sont q\ et q 2 , les équations de Lagrange s’écrivent : 

d_ r oc 

dt dqi 


OC dV _ F 
dq\ dq\ 91 


d f OC 


dt [dq 2 

Dans cette expression F* et F q2 sont les forces généralisées conjuguées des coordonnées 
généralisées respectives q\ et q 2 . Elles sont respectivement définies par 

, dans cette expression ÔW\ représente le travail des forces extérieures 

<5<?i^0 
<5<j2=0 

pour une variation 8q\ de la coordonnée qi, lorsque ôq 2 = 0. 


, dans cette expression dW 2 représente le travail des forces extérieures 

5çi=0 
<5q 2 ^0 

pour une variation ôq 2 de la coordonnée q 2 , lorsque ôq± = 0. 


_ ÔW 2 
àq-2 


F 

Q1 Sq i 


dC dV _ F 
dq 2 + dq 2 92 


5.2 Système masses-ressorts-amortisseurs 



Système à deux degrés de liberté en oscillations forcées. 


Pour étudier les particularités des oscillations forcées des systèmes à deux degrés de liberté, 
étudions le système symétrique de la figure ci-dessus, soumis à une force horizontale F appliquée 
à la première masse. 
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5.2.1 Equations différentielles 

Les équations différentielles du mouvement s’écrivent : 

mx i + (k + K) x'i + ax\ — Kx 2 = F 


—Kx 1 + mi2 + (A: + IL) X2 + 0x2 = 0 

5.2.2 Etude du régime permanent sinusoïdal 
Solution permanente 

La solution générale de système d’équations différentielles est égale à somme de la solution du 
système homogène et d’une solution particulière. La solution de l’équation homogène, en raison 
de l’amortissement, tend vers zéro lorsque le temps augmente. Lorsque le régime permanent 
s’établit, la solution devient égale à la solution permanente et s’écrit alors : 

x± = X\ cos (fit + <f> 1) 


x 2 = X 2 cos (fit + 4 > 2 ) 

Pour calculer les amplitudes X\ et X 2 , ainsi que les phases 4>i et (f> 2 , utilisons la méthode des 
nombres complexes. On peut ainsi écrire : 

x\ = Re (Xy e int ) x 2 = Re (X 2 e mt ) F = Re (F é m ) 

Dans ces expressions les amplitudes complexes sont définies par 

X { = Xi e i4>1 X 2 = X 2 e^ 2 F = F 0 e i0 

Dans ce cas les équations différentielles se transforment en équations algébriques : 

f (k + K - mfi 2 + iaiï) X x - KX 2 = F 

l -KX 1 + (k + K - mfl 2 + iaü) X 2 = 0 

Amortissement négligeable 

Considérons d’abord le cas d’un amortissement suffisamment faible pour que l’on puisse 
considérer que a ~ 0. Le système d’équations différentielles s’écrit alors 

f ( k + K- mfi 2 ) X-y - KX 2 = F 

l -KXy + (k + K - mÜ 2 ) X 2 = 0 

Les solutions de ce système sont : 

_ f (n 2 A - n 2 ) 

1 m (uf — D 2 ) (oj 2 — D 2 ) 

Y _ K F 1 

2 m 2 {ojy - D 2 ) (u 2 - D 2 ) 
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T , . k k+2K , , . 

Les pulsations u 1 = \ — et oj 2 = \ sont les pulsations propres calculées au chapitre 

V rn V rri 

précédent. La valeur de la pulsation fl^ est : 


Ü.A = 


k + K 


rn 


Les amplitudes des déplacements X\ et X 2 sont alors données par 

F 0 \n 2 A -n 2 \ 


X! = 


m 


uj 2 ~n 2 \\uj 2 -n 2 \ 


x 2 = 


K F) 


m 2 uj 2 — fl 2 | |u;| — F, 2 


Les variations des amplitudes X] et X 2 sont représentées sur les figures ci-dessous. 




Variation de X\ en fonction de fl Variation de X 2 en fonction de F, 

On remarque que le phénomène de résonance se produit pour Xi comme pour X 2 lorsque la 
pulsation d’excitation 17 est égale à l’une des pulsations propres u\ ou ui 2 du système. L’amortis- 
sement étant très faible, les amplitudes à la résonance sont très importantes. Lorsque la pulsation 
fl devient très grande, ces amplitudes tendent vers zéro. Enfin lorsque fl = l’amplitude X\ 
est égale à zéro ; pour cette raison, la pulsation fl A est appelée pulsation d’antirésonance. 


5.3 Impédance 

Considérons le système à deux degrés de liberté étudié dans le paragraphe précédent dans 
lequel nous supposons que l’amortissement est nul (a ~ 0). En régime stationnaire, on obtient 
pour l’amplitude complexe de la vitesse X_ 1 : 

+ ,-g r 

- 1 m(ï!2_„2) (!î2- w 2)- 

On en déduit l’impédance d’entrée : 

_ F _ .m (fl 2 - uj 2 ) (fl 2 - cul) 

~ E ~ ~ 1 f7 fl 2 — ü 2 a 

Les figures ci-dessous donnent les variation de X et Ze en fonction de fl. On note le phé- 
nomène de résonance lorsque la pulsation d’excitation fl est égale à l’une des deux pulsations 
propres uj\ ou co 2 . A ces pulsations, le module de l’impédance d’entrée est nul. Enfin, lorsque fl 
est égale à la pulsation d’antirésonance fl^, la vitesse de la première masse est nulle et le module 
de l’impédance d’entrée est infini. Lorsque fl — > 00 , Ze — mfl. 
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5.4 Application 

Le phénomène d’antirésonance peut être avantageusement utilisé pour supprimer une vibra- 
tion résultant d’une résonance dans un système mécanique. 



Etouffeur de vibrations . Variation de Xi en fonction de O. 


Considérons le système à deux degrés de liberté de la figure ci-dessus. Les équations différen- 
tielles du mouvement s’écrivent 

rnx i + axi + ( k + K) x\ — Kx 2 = F 

—Kx 1 + mx 2 4- Kx 2 = 0 
En régime permanent sinusoïdal, on obtient 

F 0 


Ki = 


O 2 _ K 
“ M 


X 2 = - 


m [-04 + 02 (M + f ) - £* ] + *^0 (02 - K ) 

VF,, 1 

Mm [-04 + fi 2 (*±* + f ) - ^ (02 - K ) 

Lorsque la pulsation de la force excitatrice est égale à 0^ = la masse m est immobile 

i-V, 0). 
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Si on choisit K et M telles que ^ (c’est-à-dire telles que uq = IÏa), la masse m est 

immobile lorsque la pulsation excitatrice Q est égale àwo = \ = \J jg. Dans ces conditions, 

l’ajout de M et K permet d’annuler la vibration de m à cette pulsation. Un tel dispositif constitue 
un "étouffeur” dynamique de vibrations. 


5.5 Exercices 


Exercice 1 : Dans la figure ci-dessous, M et R représentent respectivement la masse et le 

rayon de la poulie, x\ et X 2 les écarts des deux masses par rapport à leur position d’équilibre. On 
prend : M = 2(rri2 — m\) avec rri 2 = m et ko = k\ = &2 = k. Le système sera étudié en régime 
permanent avec : F(t) = Fq cos (ut). 



1. Calculer l’impédance d’entrée. 

2. Calculer les vitesses x\ (t) et X 2 (t)- 

3. (a) A quelle pulsation la masse m\ reste-t-elle immobile? 

(b) Dans ce cas, quelle est l’amplitude d’oscillation de la masse rri 2 ? 

4. (a) A quelles pulsations, la vitesse de la masse m-i est en phase avec la force F(t) ? 

(b) En comparant ces pulsations à celles calculées dans l’exercice 4 du chapitre précédent, 
déduire le déphasage de X 2 par rapport à F(t) pour chacune de ces pulsations. 

Exercice 2 : 

1. Etablir les équations différentielles régissant le fonctionnement des systèmes représentés 
par les figures 1 et 2. 

2. Pour 

- (jü = ^hk/Ami + g /r et m\ = m .2 pour la figure 1 
-et uj = \fk 2 jrn. avec m = M/2 + rri 2 pour la figure 2, 
calculer en régime permanent sinusoïdal : 

(a) l’impédance d’entrée, 

(b) 2 /i , ÿ 2 , yi et 2/2 • 
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Exercice 3 : Pour les systèmes mécaniques des figures 1 et 2 et avec F(t) = Fq cos(c ot) et 
pour lo = \Jk\/m\ , calculer : 

1. L’impédance d’entrée , 

2. La puissance instantanée fournie par le générateur mécanique et la puissance instantanée 
dissipée par le système. Comparer et commenter. 

3. La puissance moyenne fournie par le générateur mécanique et la puissance moyenne dissipée 
par le système. Comparer et commenter. 



figure 1 figure 2 
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Chapitre 6 


Généralités sur les phénomènes de 
propagation 

6.1 Propagation à une dimension 

6.1.1 Equation de propagation 

Dans les phénomènes vibratoires traités dans les chapitres précédents, nous nous sommes 
intéressés à des phénomènes ou des grandeurs physiques qui dépendaient d’une seule variable, le 
temps. Nous allons maintenant examiner toute une série de phénomènes qui sont décrits par une 
fonction qui dépend à la fois du temps t et d’une variable d’espace, x par exemple. 

Ces phénomènes sont régis par une équation aux dérivées partielles, appelée équation de 
d’Alembert ou équation d’onde ou encore équation de propagation à une dimension de la forme : 

d 2 s 1 d 2 s 

dx 2 V 2 dt 2 

dans laquelle V est une grandeur physique qui a les dimensions d’une vitesse et sera appelée dans 
la suite vitesse de propagation. 

6.1.2 Solution de l’équation de propagation 
Méthode de d’Alembert 

Pour résoudre l’équation des ondes à une dimension, opérons le changement de variable 
suivant : 


7/ = t — 


X 

V 


Ç 


t + 


x 

V 


Calculons les dérivées partielles par rapport à t et x, en fonction des dérivées partielles par 
rapport à rj et £. 

Sachant que : 

drj 

dt^dt^ 

et que 

dr) dÇ 1 

dx dx V 

on obtient 
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d s ds di i d s <9£ ds ds 

dt di] dt. dt i dt di] <9£ 


ds 

ds drj ds <9£ 1 

ds 

ds~\ 

dx 

di] dx <9£ dx V 

di] 

I 


En tenant compte de ces résultats et sachant que 


d 2 s d 2 s 

dr]dt; dÇdr] 


on obtient : 

d 2 s d 2 -s 9 d 2 s d 2 s 

dt 2 dif ^ dr/dÇ <9 £ 2 

d 2 s 1 d 2 -s ^ d 2 s d 2 s 

dx 2 V 2 dr/ 2 ^ dr/dÇ dÇ 2 

En remplaçant dans l’équation d’onde et par les expressions ci-dessus, on obtient l’équa- 
tion d’onde exprimée en fonction des dérivées partielles par rapport aux variables i] et £ : 

d 2 s 

drjdÇ 


Cette dernière équation peut s’écrire 

d ds 
dt i dt] 

Un intégration par rapport à £ donne : 


où / ( ï ] ) est une fonction qui ne dépend que de rj (et pas de £). Enfin une intégration par rapport 
à r) donne : 

s(r,,0 = F( V ) + G(0 

où F (rj) , qui ne dépend que de r], est une primitive de / ( 77 ). La fonction G (£) est une fonction 
qui ne dépend que de £. En revenant aux variables x et t, on obtient la solution générale de 
l’équation des ondes à une dimension : 

s(x,t) = F (t- +G (t+ 

Les fonctions F [t — y") et G (t + y) sont des fonctions dont la nature est fixée par les conditions 
aux frontières imposées à la solution s (. x,t ). 
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Propriétés des solutions particulières F (t - f ) et G (t + f ) 

Propriétés de F (t — |)) On étudie le cas de la solution particulière F (t — y). Pour cela 
on suppose que les conditions aux frontières sont telles que G (t + y) est constamment nulle. 
On considère à l’instant t\ un point d’abscisse x\ . La valeur de la fonction s en ce point et à cet 
instant est s(aq,ti). On recherche à un instant t 2 postérieur à t\ ( £2 > fi) la position X2 d’un 
point pour lequel la valeur de s est la même que la valeur qu’elle avait en x\ à l’instant t\. Ce 
problème est formulé par l’égalité suivante : 

S (.Tl, fl) = S (x 2 ,t 2 ) 


Ce qui se traduit par 


Cette équation est satisfaite si 


Tl 


V 


F (h — rrf) = F (t 2 - ~=r 


X 2 


V 


h - 


Tl 

V 


h ~ 


T2 

V 


D’où la valeur de x 2 : 


X 2 = Tl + V (' t 2 ~ 1 1 ) 

Comme t 2 > ti, x 2 est supérieure à ti et ces deux points sont distants de 


x 2 — Tl = V (t 2 - 1 1) 

F {t ~ y) correspond à une onde se propageant dans le sens des x croissants (Voir la figure 
ci-dessous). F (t — y) est appelée onde progressive et cette expression constituera dans la suite 
la définition d’une onde progressive. 



Direction de 
propagation 




t=t 2 >ti 



xa-x^Vfe-t]) 

Onde progressive dans le sens des x croissants : F (t — y) 


Propriétés de G (t + y ) On étudie le cas de la solution particulière G (t + y). Pour cela 
on suppose que les conditions aux frontières sont telles F (t — y) est constamment nulle. On 
considère à l’instant t\ un point d’abscisse ti . La valeur de la fonction s en ce point et à cet 
instant est s(Ti,ti). On recherche à un instant t 2 postérieur à t\ (t 2 > ti) la position x 2 d’un 
point pour lequel la valeur de s est la même que la valeur en x\ à l’instant t,\. Ce problème est 
formulé par l’égalité suivante : 

s (Tl,fl) = -s {x 2 ,t 2 ) 

Ce qui se traduit par 
Cette équation est satisfaite si 


G{tl + v) - G { t2 + V 


Tl T 2 

h + y=i 2 + y 
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D’où la valeur de x 2 ■ 


x 2 = xi - V (t 2 - h) 


Comme t 2 > t\, x 2 est inférieure à x\. Ces deux points sont distants de 


xi - x 2 = V (t 2 - ti) 

G (t + y) correspond à une onde se propageant dans le sens des x décroissants (Voir la figure 
ci-dessous). 


c 


Direction de 
propagation 



t=t 2 >ti 



Xi-x 2 =V(t 2 -t,) 


Onde progressive dans le sens des x décroissants : G [t + y) 


6.1.3 Onde progressive sinusoïdale 

On considère une onde progressive se propageant dans la direction de l’axe des x, telle que 
le point d’abscisse x = 0 est soumis à une vibration sinusoïdale de la forme 


s (x = 0, t) = S 0 cos (ut) 


Le point se trouvant à l’abscisse x > 0 aura la même vibration que celle du point x = 0 mais 
avec un retard égal à y : 


s ( x , t ) = S 0 cos 


, x 

u[t -y 


Cette expression constitue la définition d’une onde progressive sinusoïdale (ou harmonique) ; elle 
peut être écrite sous la forme : 


s ( x , t ) = Sq cos [ut — 4> (^)] 


où 4> (x) = yx représente le déphasage lié au temps de propagation y. On dit que cj) (x) représente 
le déphasage dû à la propagation. L’onde progressive sinusoïdale s’écrit sous la forme suivante 
qui permet de mettre en évidence la double périodicité (dans le temps et dans l’espace) : 


s (x, t ) 


Sç> cos 




La quantité T = ^ est la période temporelle tandis que la quantité À = VT est la longueur 
d’onde qui constitue la période spatiale. On peut vérifier aisément que : 


s(x,t + nT) = s(x,t ) 
s (x + nX, t) = s(x,t ) 


où n est un nombre entier. 
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L’onde progressive s’écrit souvent : 

s (. x , t ) = S 0 cos [ut — kx\ 

où k = y = ^ est appelé le module du vecteur d’onde qui s’exprime en mT 1 . 

On utilise très souvent la notation complexe d’une onde progressive sinusoïdale : 

s(x,t) = S 0 
s (x, t) = S e iujt 

où 5 = So e~ lkx représente l’amplitude complexe de l’onde progressive sinusoïdale. Le module 
So de S est l’amplitude de l’onde tandis que son argument — kx représente le déphasage dû à la 
propagation. 

6.1.4 Superposition de deux ondes progressives sinusoïdales 

Cas de deux ondes de même fréquence se propageant dans le même sens 

Considérons deux ondes de même fréquence et de même direction de propagation, d’ampli- 
tudes respectives Si et S 2 -, et de phases respectives 4>\ et <f> 2 - L’onde résultante sera alors : 

g f . x ; _ J 5 1 e i{ojt-kx+<jn) ^ e i(ujt-kx+<f> 2 ) _ g & i(wt-kx+<j>) 

ou encore en notation réelle : 


avec 


s ( x , t) = S cos (ut — kx + 4>) 


S= y / Si + + 25i5 2 cos (</>! - <fo) 

et 


= Arct ( Si sin(0i) + 5 2 sin (0 2 ) \ 

‘ 9 \S! cos (</>!) + S 2 cos(cf>2)J 

La superposition de deux ondes harmoniques de même fréquence, et qui se propagent dans la 
même direction, donne une autre onde harmonique progressive de même fréquence, d’amplitude 
S et de phase <j>. 


Cas de deux ondes de même fréquence se propageant dans des sens opposés 

Si par contre, on superpose deux ondes harmoniques de même fréquence mais se propageant 
dans des sens opposés, le résultat est tout autre. En effet, dans ce cas : 


s (x, t ) = Si + S 2 e^+k'+fc) 


Sié^ 1 e~ ikx + S 2 e^ 2 e +ikx 


e iu}t 


et on ne peut plus écrire l’onde résultante sous la forme d’une onde progressive simple. Un cas 
particulier important se produit quand les deux amplitudes sont identiques. Si on note : 


Si = S 2 = So 


on a : 


s (x, t ) = 2Sb cos 



4 > 1 - 4>2 
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et donc en notation réelle : 

s (x, t ) = 2So cos ^ kx + c h__Él' S j cos _)_ < ^ >1 ^ ^ 

Ce mode de vibration est très différent d’une onde progressive puisque tous les points x de la 
corde vibrent en phase avec des amplitudes différentes. En particulier, il existe une série de 
points : 


avec 



< fil - 4 > 2 

2tt 


A 

2 


n = 0, ±1, ±2, • • • 

où l’amplitude de vibration est constamment nulle. On dit dans ce cas que l’onde est sta- 
tionnaire et que les points x n sont les nœuds de l’onde. Entre chaque paire de nœuds existe un 
ventre où l’amplitude de vibration est maximum et égale à 2Sq. On note aussi que l’intervalle 
entre deux nœuds est égal à une demi-longueur d’onde A/2. 


6.1.5 Vitesse de phase 


Considérons une onde progressive sinusoidale qui se propage dans le sens des x croissants. 
Un point d’abscisse x possède, à l’instant t, l’élongation : 

s ( x , t) = S 0 cos (ut — kx) 

Entre l’instant t et t + At l’onde progresse d’une quantité Ax. A l’instant t + A t, le point 
d’abscisse x + Ax possède la même élongation que celle que possédait le point d’abscisse x à 
l’instant antérieur t. Ceci se traduit par l’égalité : 

s ( x , t) = s (x + Ax, t + At) 

S 0 cos (ut — kx) = Sq cos [u ( t + At) — k (x + Ax)] 

Cette égalité est satisfaite si les phases instantanées sont égales : 

ut — kx = u (t + At) — k (x + Ax) 


Soit encore 

On définit la vitesse de phase V/, 


u At = k Ax 

qui s’exprime en fonction de te et A: par : 



Si la vitesse de phase ne dépend pas de te, le milieu est dit non dispersif. Dans le cas contraire il 
est dit dispersif. 

La figure ci-dessous permet d’illustrer la notion de vitesse de phase en considérant deux 
représentations à des instants différents d’une corde parcourue par une onde . La courbe continue 
représente l’ensemble des points de la corde à l’instant t. Le point de la corde d’abscisse x est 
représenté par le point blanc, tandis que le point d’abscisse x + Ax est représenté par le point 
noir. On constate qu’entre les instants t et t + At chacun de ces points suit une trajectoire 
rectiligne et le déplacement du point noir à l’instant t + At est égal au déplacement du point 
blanc à l’instant t. En particulier la crête de la corde, correspondant à une valeur particulière 
de la phase instantanée, semble se déplacer dans le sens de propagation de l’onde avec la vitesse 
de Vfj, mais la trajectoire de chaque point matériel est une trajectoire rectiligne perpendiculaire 
à la direction de propagation. 
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Onde progressive sinusoïdale 

Trajectoire de deux points o et • entre les instants t et t + At 

6.1.6 Vitesse de groupe 

La vitesse de phase n’est pas nécessairement la vitesse que l’on observe lorsqu’on analyse 
un mouvement ondulatoire. En général une onde n’est pas parfaitement sinusoïdale mais a une 
durée limitée et se présente sous la forme d’un train d’onde appelé communément ”pulse” ou 
"groupe” qui se propage avec une vitesse Vq appelée vitesse de groupe. Cette onde sous la forme 
d’un puise contient plusieurs fréquences. Si la vitesse de phase est indépendante de la fréquence 
(Milieu non dispersif) alors toutes les fréquences qui constituent le puise se propagent à la même 
vitesse et le puise se propage avec une vitesse de groupe égale à la vitesse de phase. Mais si le 
milieu est dispersif (i.e. la vitesse de phase dépend de la fréquence), alors le puise se propage 
avec une vitesse de groupe différente de la vitesse de phase. 

Pour illustrer ce phénomène, considérons une onde constituée de deux ondes de fréquences 
différentes et de même amplitude. En x = 0, cette onde s’écrit par exemple sous la forme : 

S (0, t ) = Sq COS (u\t) + Sq COS (u> 2 t) 

Cette onde peut s’écrire encore : 

s (0, t ) = 2Sq cos (c ost) cos (ut) 


où 

UJ2 — Ui ; U 2 + Ui 

u, B = — et u = — 5 — 

Si u\ est voisine de U 2 , la vibration résultante se présente sous la forme d’une sinusoïde de 
pulsation u dont l’amplitude est modulée par un battement de pulsation ub (Modulation d’am- 
plitude). 

En un point x > 0, l’onde obtenue résulte de la superposition de ces deux ondes qui se sont 
propagées à des vitesses différentes car le milieu de propagation est supposé dispersif : 

s ( x , t) = Sq cos (a )\t — k\x ) + Sq cos (a;2 1 — k2x) 


s (x, t ) peut s’écrire : 


s (x, t) = 2 Sq cos (ust — ksx) cos (ut — kx ) 


Dans cette expression : 


kB 


k2 - h 
2 


et k 


k 2 + ki 
2 
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L’amplitude du battement se propage à une vitesse qui est la vitesse de groupe définie par la 
relation : 


Vg = 


ub 
k b 


ÜÜ2 ~ Wi 

k 2 - h 


Comme k 2 est peu différent de k \ . la vitesse de groupe est définie par : 


V G 


lim 

k 2 — >ki 


0J 2 ~ lü 1 

k 2 - ki 


du; 

dk 


Tandis que la sinusoïde contenue à l’intérieur du battement se propage à la vitesse de phase : 




Propagation d’un paquet d’ondes : Les flèches verticales noires correspondent au maximum des 
battements qui se propagent à la vitesse de groupe. Les flèches verticales blanches 
correspondent au maximum des vibrations qui se propagent à la vitesse de phase. 

6.1.7 Onde vectorielle 

Dans ce qui précède, la quantité s (. x , t) représente une grandeur scalaire, mais certains phé- 
nomènes décrits par des vecteurs conduisent à des équations similaires : 

d 2 Â 1 d 2 Â _ 

dx 2 V 2 d t 2 

Le vecteur A. défini dans un milieu à trois dimensions, possède trois composantes A x , A y , A z et 
l’expression ci-dessus signifie que chacune de ces composantes satisfait individuellement l’équa- 
tion de propagation : 




Chacune de ces composantes A x , A y , A z se propage en t — y et t + y. 


d 2 A r 

1 

d 2 A r 

= 0 

dx' 2 


dt' 2 

d 2 Ay 

1 

dAy 

= 0 

dx 2 


dt 2 

d 2 A z 
dx 2 

i 

V 2 

d 2 A z 
dt 2 

= 0 
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6.2 Propagation dans l’espace à trois dimensions 

6.2.1 Equation de propagation 

Dans un système de coordonnées cartésiennes, l’équation de propagation dans l’espace à trois 
dimensions s’écrit sous la forme : 

d 2 -s d 2 -s d 2 s 1 Ô 2 s 

dx 2 dy 2 dx 2 V 2 dt 2 

On définit le laplacien scalaire de s par l’expression ci-dessous : 

d 2 s d 2 s d 2 s 

dx 2 dy 2 dx 2 


et l’équation des ondes s’écrit sous la forme condensée : 


As 


1 d 2 s 

v^dt 2 = 


s est fonction du temps mais également des cordonnées du point M où la fonction 
calculée. Si l’on appelle le vecteur position r = OM, la quantité s dépend du temps et 
position r ; on écrit s ( r,t ). 


s doit être 
du vecteur 


6.2.2 Onde plane progressive sinusoïdale 
Définition 

L’onde progressive sinusoïdale (ou harmonique), se propageant dans une direction donnée 
par un vecteur unitaire u est définie par : 

s (r, t) = S 0 cos (ut 

où le vecteur k = k ü est appelé le vecteur d’onde. Si les composantes du vecteur k sont k x , k y 
et k z , alors l’onde plane est définie par 

s (x, y, z , t ) = So cos (ut — k x x — k y y — k z z) 

On peut utiliser la notation complexe pour représenter l’onde plane progressive sinusoïdale qui 
s’écrit dans ce cas sous la forme : 


s(f,t) = S 0 e < ( wt -*' ? ) 


Relation de dispersion 


En remplaçant s(r, t ) par son expression dans l’équation de propagation, on obtient la relation 
k = k (u) qui est la condition pour que l’onde plane définie ci-dessus constitue une solution 
particulière de l’équation d’onde. Cette relation est appelée la relation de dispersion et elle 
s’écrit : 


k = 


u 

V 
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Surface d’onde 

On appelle surface d’onde ou surface équiphase, l’ensemble des points de l’espace pour les- 
quels, au même instant, s (r, t) a la même valeur. Recherchons la surface d’onde passant par un 
point Mo à un instant t ; cette surface est l’ensemble des points M de l’espace pour lesquels 
l’ égalité suivante est satisfaite : 

s (r, t) = s (r 0 , t ) 

Cette égalité se traduit par : 

S 0 = s 0 e i(w*-fc-ro) 

Cette égalité est satisfaite si 

k ■ f = k ■ ro 

Tenant compte des propriétés du produit scalaire, on obtient 


k ■ (r — ro) = 0 


La surface d’onde passant, à l’instant t , par le point Mo est l’ensemble des points M satisfaisant 
l’équation ci-dessus. Cette surface est un plan passant par Mo, et perpendiculaire à la direction 
du vecteur d’onde, donc à la direction de propagation. L’onde est dite plane. 



Il existe d’autres types d’ondes définis par les surfaces d’onde respectives : par exemple les 
ondes sphériques pour lesquelles les surfaces d’onde sont des sphères ou les ondes cylindriques 
pour lesquelles les surfaces d’onde sont des cylindres. 

Polarisation 

Dans le cas d’une onde plane progressive sinusoïdale représentée par une quantité vectorielle 
A(r,t), cette quantité peut avoir différentes orientations par rapport aux surfaces d’ondes : 

1. A est constamment perpendiculaire à la surface d’onde, ou de manière équivalente parallèle 
à la direction de propagation : l’onde est dite longitudinale. 

2. A est contenu dans la surface d’onde, ou de manière équivalente perpendiculaire à la di- 
rection de propagation : l’onde est dite transversale. Dans ce cas, l’extrémité du champ 
vectoriel A peut décrire une trajectoire rectiligne : l’onde transversale est dite à polarisa- 
tion rectiligne. Elle peut décrire une trajectoire circulaire (onde transversale à polarisation 
circulaire), ou une trajectoire elliptique (onde transversale à polarisation elliptique). 
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6.3 Exercices 

Exercice 1 : Vérifier que les fonctions suivantes : 

1. u (x,t) = Asm [u [t — y)) 

2. u (x, t) = A cos (k {x — Vt)) 

3. u (x, t ) = a(x + Vt) 2 

4. u ( x , t) = A exp ( iuj [t — y) ) 

5. u (. x , t) = A exp ( iuj (t + y) ) 

6. u ( x , t) = A exp (ica (t — y') ) + B exp (ica (t + y)) 

sont solutions de l’équation : 

d 2 u 1 d 2 u 

dx 2 V 2 dt 2 

où x, t et V représentent respectivement la position, le temps et la vitesse de propagation. 
Déterminer les dimensions des constantes A, ca, k et a. 

Exercice 2 : On étudie la propagation d’un ébranlement transversal sur une corde. Cet ébran- 
lement se propage dans le sens des x croissants avec une vitesse V. A l’instant t = 0 s, la forme 
de la corde est donnée par : u(x, 0) = Ae~ ax . Donner l’expression de la forme u{x,t) de la 
corde à un instant t. Représenter graphiquement cette corde aux instants t = 0 s, et t = 0.3 s 
pour : A = 1 cm, a = 0.5 cm -2 , V = 20 cm. s -1 . 

Exercice 3 : Le dispositif représenté sur la figure ci-dessous permet de communiquer à l’ex- 
trémité S d’une corde tendue horizontalement, une vibration verticale sinusoïdale entretenue : 
u (0, t) = 5 sin (c ot) (en cm) 



1. La longueur totale de la corde est L = 5 m et sa masse est m = 100 g. Sachant que la 
vitesse de propagation des ondes transversales dans une corde de masse linéique fi et tendue 
par une tension T est donnée par : V = \/T / fi , quelle doit être la valeur de la masse M 
qui tend la corde pour que la longueur d’onde des vibrations transversales de pulsation soit 
A = 1 m lorsque le moteur tourne à la vitesse de 3600 tr/mn. 

2. Juste avant la poulie, la corde est pressée à frottement doux entre deux plaques de feutre, 
empêchant toute réflexion de se produire. 

(a) Ecrire l’expression en fonction du temps de l’élongation u(x,t ) du point A tel que 
SA = x. 

(b) Donner les abscisses des points qui vibrent en phase avec S et de ceux qui vibrent en 
opposition de phase avec S. 

(c) Représenter sur un même graphique le mouvement de S entre les instants t = 0 et 
t = 1/30 s et le mouvement du point A tel que x = 2, 75 m. 

(d) Représenter l’aspect de la corde sur un même graphique aux instants t = 0 et t = 
1/600 s, entre x = 0 et x = 4 m. 
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Chapitre 7 


Cordes vibrantes 


7.1 Equation des ondes 

Considérons une corde tendue, rectiligne selon la coordonnée x. et de longueur infinie. Nous 
allons étudier la propagation d’un faible ébranlement le long de la corde. Supposons que cet 
ébranlement se produise suivant l’axe 0 y. 

Etudions l’équation du mouvement de cette corde. Nous dénoterons par T la tension à laquelle 
est soumise la corde. On considère en un point d’abscisse x un segment très court de cette corde, 
de longueur Ax. La masse Am du segment est donnée par : 

Am = /iAx 

où /i est la densité linéique de masse de la corde, c’est-à-dire la masse par unité de longueur 
qui s’exprime en kg/m. 



Corde vibrant transversalement. 


Dans une situation hors équilibre, le segment n’est plus droit, il présente une courbure. Nous 
considérons des mouvements d’oscillation de la corde de petite amplitude 


u (x, t) = u ( x , t ) e y 

si bien que nous pouvons faire l’approximation : 

3v 

sin Wl* = tan WU = fa 

du 

sin (0)UAz = tan ( 0 )lz+Ax = fa 


x-\-Ax 
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Cette approximation néglige aussi l’allongement du segment, et considère donc la tension T 
comme constante. La force appliquée sur le segment dans la direction y est la résultante de la 
force appliquée au point x, qui est une force appliquée vers le bas et égale en module à 

Ou 

F (x, t) = T sin (6)\ X = T tan (6)\ X = T — 
et de la force appliquée au point x+Ax qui est vers le haut et égale à 


F (x + Ax, t ) 


T sin ( d )\x+Ax = T tan ( 0 )lx+A* = T 


du 

dx 


x-\-Ax 


La force totale dans la direction y est donc : 


R = F (x + Ax, t ) 


F fat) = T 


du 

du 


dx 

x+Ax dx 

X- 


rn d 2 u 

T 

dx 2 


Ax 


Nous pouvons appliquer maintenant la loi fondamentale de la dynamique au segment Ax. La 
force dans la direction y doit être égale au produit de la masse Am du segment par l’accélération 
de celui-ci. Donc : 


R = A m 


d 2 u 
dt 2 


d 2 u d 2 u 

^ dt 2 dx 2 


Si on définit V 



qui a la 


dimension d’une vitesse, on constate que : 


d 2 u 1 d 2 u 

dx 2 V 2 dt 2 

qui est l’équation d’onde de la corde. V est la vitesse de propagation de cette onde. 


7.2 Ondes progressives harmoniques 

7.2.1 Définition 

Une onde progressive harmonique se propageant selon Ox est définie par : 

u (x, t ) = Uo cos (ut — kx ) 
ou encore en notation complexe 


u fat) = U 0 e i(ujt ~ kx) 

où k = y = ^ est le module du vecteur d’onde, À étant la longueur d’onde. 

7.2.2 Force en un point 

On appelle force en un point, la projection selon Oy de la force exercée, en ce point, par la 
partie gauche de la corde sur la partie droite : 
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Dans le cas d’une onde progressive sinusoïdale, cette relation devient : 

F(x, t ) = ikT U 0 e *( art ~ fcx ) 

La vitesse de particules s’écrit : 

ù(x,t) = ^=icjU 0 

On constate que pour une onde progressive la vitesse de particules ù est en phase avec la 
force F. 

7.2.3 Impédance 

On appelle impédance en un point le rapport de l’amplitude complexe de la force à l’amplitude 
complexe de la vitesse de particule 

Z{X)= Ù L 

Dans le cas d’une onde progressive, on obtient : 

Z (x) = gV = T 

La quantité y/jZT définit l’impédance caractéristique de la corde 

Z c = ^/JZr = gV 

On obtient une propriété de l’onde progressive plane 

Z (x) = Z c Vx 

7.3 Oscillations libres d’une corde de longueur finie 



Corde de longueur L fixée aux extrémités 


Considérons une corde de longueur L fixe aux points x = 0 et x = L. Recherchons une 
solution de l’équation d’onde sous la forme : 

u (x, t) = g (x) / (t) 

En remplaçant dans l’équation de propagation, on obtient : 

1 d 2 g _ 1 1 d 2 f 
g dx 2 V 2 f dt 2 

Le membre de gauche de cette équation ne dépend que de x, tandis que le membre de droite 
ne dépend que de t. Ces deux expressions sont donc égales à une constante qui doit être un 
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nombre réel négatif que nous posons égal à —k 2 car la solution ne doit pas tendre vers l’infini 
lorsque t tend vers l’infini. Posons u = k V . On en déduit que : 


d 2 g 

dx 2 


= ~k 2 g 




dt 2 


Les solutions de ces deux équations différentielles sont de la 


fnr m p 


/ = A cos (ut) + B sin (ut) 


g = C cos (kx) + D sin (kx) 

La solution de l’équation d’onde peut alors s’écrire sous la forme : 


u (x, t ) = [A cos (ut) + B sin (ut)] [C cos (kx) + D sin (kx)] 
Tenant compte des conditions aux limites 

u ( 0 , t) = 0 


u (L, t) = 0 

on obtient 


C = 0 

k = n J où n = 0, 1, 2, • • • 

La solution de l’équation d’onde qui satisfait ces conditions aux limites est donc la somme 
d’une infinité de termes : 


avec 


OO 

u (x, t) = ^2 [a n cos (u n t) + b n sin (u n t)] sin (k n x) 

n = 0 


7 r 


k n = n— et u n 

lu 


k n V = n 


nV 

X 


Les u n sont les pulsations propres. Les coefficients a n et b n sont déterminés par les conditions 
initiales du mouvement. Supposons qu’à t = 0 nous imposions à la corde une certaine forme 
initiale u(x, 0) = uq(x) et une vitesse initiale 


il (0, x) = vq (x) 

Dans ce cas nous aurons les conditions initiales suivantes : 


OO 

u o (x) = J2 a n sin (k n x) 

n = 0 


oo 

v 0 (x) = -Unbn sin (knx) 
n = 0 

On doit inverser ces équations pour obtenir les coefficients a n et b n . La méthode de Fourier 
consiste à les multiplier par sin (k m x) et les intégrer entre 0 et L. Si on utilise les intégrales : 


{ 0 si m 7^ n 

L 

— si m = n 
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on obtient 


a n = — f uq (x) sin ( n — dx 
LJ o \ L / 

J vq ( x ) sin ^n— ^ d‘ x 


hn — 


'n 1 - 1 J 0 


7.4 Réflexion et transmission 

7.4.1 Réflexion et transmission entre deux cordes semi- infinies 

■=> 


incidente réfléchie 


transmise 


T, Mi 


0 


T, M 2 

r 2 -p 

I Ml U2 

Z 1 = V/¥ z 2 = 

Réflexion transmission dans deux cordes semi-infinies 

Soit deux cordes de longueur semi-infinie, reliées en x = 0. Leurs masses linéiques sont 
respectivement ji\ et /x 2 - Lorsqu’une onde venant de — oo se propage vers x = 0 dans la première 
corde, elle donne naissance au point de jonction, x = 0, à une onde réfléchie et une onde transmise. 
L’écriture de la continuité du déplacement et de la force en x = 0 permet d’obtenir le coefficient 
de réflexion R u et le coefficient de transmission T u définis respectivement par : 

„ _Vr 

r --1a 

T _Ut 

u Ut 

où Ui , Ur et Ut sont les amplitudes des déplacements associés respectivement à l’onde 
incidente, l’onde réfléchie et l’onde transmise. On en déduit : 


Rv = 


Z\ — Z-2 
Z\ + Z-2 


T u = 


2Z 1 
Z\ + Z 2 


7.4.2 Réflexion sur une impédance quelconque 



Corde semi-infinie terminée par une impédance Z t 
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Soit une corde de longueur semi-infinie, de masse linéique //. tendue horizontalement avec une 
tension T et terminée en x = 0 par une impédance mécanique Zt ■ Lorsqu’une onde harmonique 
se propage dans la corde de — oo vers x = 0, elle subit une réflexion en ce point. Sachant que le 
déplacement de particules s’écrit : 

u (x, t) = Ui e i 2 3 ( ut ~ kx ) + U R e i ( wt+kx ) 

on en déduit la vitesse de particules et la force en un point d’abscisse x 
û = § ± = iu[Ui e i ( ut ~ kx ') + U R e *(«t+*»)] 


F = = —ikT \Ui e i ( ut ~ kx ') 


u r e ^ ut+kx y 


En x = 0, les conditions aux limites s’écrivent : 


F(0,t) 

ù(0, t ) 


On en déduit le coefficient de réflexion R u en fonction de l’impédance caractéristique Z c et 
de l’impédance Zt placée à l’extrémité de la corde : 


_ Ur _ Z c - Zr 
U i Z c + Zt 


7.5 Exercices 


Exercice 1 : Soit une corde de longueur L = 2 m, de masse totale m = 80 g et soumise 

à une tension T réglable. Comme indiqué sur la figure ci-dessous, elle est fixée en A. à un 
ressort vertical de raideur K dont l’extrémité C subit un déplacement forcé vertical sinusoidal 
de fréquence / = 50 Hz, d’amplitude so = 1 cm et décrit par : s(t ) = s^e 

La seconde extrémité B{x = L ) de la corde est fixée à un amortisseur de coefficient de 
frottement visqueux a égal à 0,2 N. s/m. Des guides parfaitement glissants n’autorisent que 
les mouvements verticaux (transverses) des points A, B et C. Le mouvement d’un point M 
quelconque d’abscisse x sur la corde est alors représenté par son élongation u y (x,t ) . La tension 
T est réglée de telle sorte qu’il n’y ait pas d’onde réfléchie en B. 


y 



1. Que représente l’amortisseur par rapport à la corde? Calculer l’impédance en un point 
quelconque de la corde. En déduire la tension T de la corde. 

2. En déduire la vitesse de phase V et la longueur d’onde . 

3. Montrer que les extrémités H et H vibrent en phase. 
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4. Quelle est la valeur de l’impédance de la corde au point A. Calculer l’amplitude et la phase 
(p de l’élongation u y (x,t ) du point A. En déduire l’expression de l’élongation u y (x,t ) d’un 
point quelconque de la corde. Dans quel cas la phase ip est nulle? Que vaut, dans ce cas, 
l’amplitude de vibration de chaque point de la corde? 

Exercice 2 : Soit deux cordes de longueur semi-infinie reliées en x = 0. La corde qui s’étend 
de — oo à 0 a une densité linéique /x i . La seconde corde qui s’étend de 0 à +oo a une densité 
linéique fi 2 = 0.25 n\. Lorsqu’une onde incidente sinusoïdale se propage de — oo dans le sens des 
x positifs, elle subit une réflexion partielle en x = 0. L’amplitude de l’onde incidente est Uq et 
sa pulsation est u . 

1. Calculer le coefficient de réflexion en x = 0. 

2. Montrer que l’onde résultante dans la première corde a une amplitude qui varie entre deux 
valeurs extrêmes U ma x et U m i n . 

(a) Donner l’expression de Umax en fonction de Uq et calculer la position des maxima en 
fonction de la longueur d’onde À . 

(b) Donner l’expression de U m in en fonction de Uq et calculer la position des minima en 
fonction de la longueur d’onde . 

(c) Calculer le taux d’ondes stationnaires. 

Exercice 3 : Une corde de longueur L est tendue entre deux points situés respectivement en 
x = 0 et x = L. Le point situé en x = 0 est fixe et le point situé en x = L est relié à un ressort 
de raideur K. La tension de la corde est T. La corde est horizontale à l’équilibre et on pourra 
négliger son poids. On étudie les ondes transverses stationnaires sinusoïdales de pulsation oj. La 
vitesse de propagation est V. 



1. Ecrire le déplacement u y (x,t) en un point quelconque d’abscisse x et à l’instant t. 

2. En écrivant la condition limite pour x = 0, montrer que u y (x,t ) peut se mettre sous la 
forme : u y (. x,t ) = Ae %wt f (x) ; où f(x) est une fonction que l’on explicitera. 

3. Ecrire la condition aux limites en x = L. 

4. Montrer que les pulsations propres doivent vérifier la condition : tg (^) = où C est 
une constante que l’on précisera. Proposer une méthode de résolution de cette équation. 

5. Déterminer les trois premières pulsations propres uj\ , W 2 et W 3 dans le cas où — > 00 . 

Exercice 4 : On considère une corde homogène de masse par unité de longueur // tendue 

horizontalement avec une tension T très grande devant le poids de la corde.. La corde de longueur 
infinie est terminée en x = 0 par un amortisseur de coefficient de frottement visqueux a . Une 
onde incidente transversale arrive de —00 et se propage dans la corde dans le sens des x croissants. 
Cette onde correspond à un déplacement transversal donné par : 

Ui(a :,t) = Ui e < ( wt - fce ) 

où k représente le module du vecteur d’onde. 
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0 




X 


□ 

□ 

a 



1. Donner l’expression du coefficient de réflexion R en amplitude de déplacement en x = 0. 
Quel est le module et l’argument de R dans le cas particulier où a < 

2. Ecrire l’expression du déplacement résultant u(x, t), en chaque point de la corde, en fonction 
des données du problème et du coefficient de réflexion R. 

3. Montrer que le déplacement résultant en chaque point de la corde peut alors s’écrire sous 
la forme de la somme d’une onde progressive et d’une onde stationnaire et qu’il s’écrit sous 
la forme : 

U (x,t) = u P +U(x ) e iut 

U p et U(x) sont deux nombres réels qui représentent respectivement l’amplitude de l’onde 
progressive et l’amplitude de l’onde stationnaire. Calculer l’amplitude U p en fonction de R 
et Ui , ainsi que l’amplitude U(x) en fonction de R, Ui , k et x. 

4. Que devient u(x, t ) dans les cas particuliers suivants : 

(a) a = y/jüT, 

(b) a = 0. 


Exercice 5 : Une corde de longueur L et de masse m est tendue entre deux points fixes avec une 
tension T. A l’instant t = 0, la corde est pincée en son milieu, écartée par rapport à l’horizontale 
d’une distance h puis lâchée sans vitesse initiale. Les positions u(x, 0) des différents points de la 
corde, à l’instant t = 0, sont représentées sur la figure ci-dessous. Etablir l’expression de u(x,t). 
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Chapitre 8 

Ondes élastiques dans les solides 


8.1 Propriétés élastiques des solides 


8.1.1 Déformation 




Allongement d’un barreau sous l’action d’une force de traction. 


Soit un solide déformable, continu et isotrope, sous la forme d’un barreau rectiligne de faible 
section, c’est-à-dire dont les dimensions latérales sont faibles devant sa longueur £, dont une 
extrémité est fixée sur un support rigide fixe. Sous l’action d’une force de traction F appliquée 
à l’extrémité libre, le fil s’allonge. La déformation persiste tant que la traction est maintenue. 
L’allongement relatif est : 

81 _£' -l 

J ~ £ 

où £' est la longueur du barreau déformé. Chaque partie du barreau n’est pas nécessairement 
déformée de la même manière. Il faut donc définir un allongement relatif local, c’est-à-dire qui 
dépend de la position en chaque point d’abscisse x. 


F (*) 

◄ 


x x + Ax 

IX 


F(x + Ax) 


\L 


u x (x) 'u x (x + Ax) 


F(x) 

◄ 

r 7 

L 

F(x + Ax) 

w. 




Déformation locale d’un barreau. 


Pour cela, considérons la portion du barreau située initialement entre les points d’abscisses 
respectives x et x + Ax. Lorsque la traction est exercée, ces points se déplacent et leurs abscisses 
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respectives deviennent x + u x (x) et x + Ax + u x (x + Ax). L’épaisseur de l’élément situé entre 
les deux sections devient : 


Ax' = [x + Ax + u x (x + Ax)] — [x + u x (x)] 


Si l’épaisseur Ax est suffisamment petite, un développement en série de Taylor au premier 
ordre permet d’écrire : 

du x 

u x (x + Ax) = u x (x) + — — Ax 

ax 

L’allongement relatif de l’élément est : 

Ax' — Ax u x (x + Ax) — u x (x) 

Ax Ax 

Par définition, la déformation e au voisinage du point d’abscisse initiale x est la limite de ce 
rapport lorsque Ax tend vers zéro, c’est-à-dire la dérivée de u x par rapport à x : 


, u x (x + Ax) — u x (x) du 

£ = l lm 7 = 77“ 

Ax— >0 Ax ox 

Il convient de noter que la déformation £ qui est le rapport de deux longueurs est sans 
dimension. Lorsque e est négative on dit que le milieu subit une contraction locale et si e est 
positive le milieu subit une extension locale. 


8.1.2 Contrainte moyenne 

Le fil s’est allongé sous l’effet de la force extérieure F exercée normalement sur sa section 
droite de surface S. Dans le cas général, cette force n’est pas constante le long du fil, elle dépend 
de l’abscisse x. On définit la contrainte moyenne par le rapport : 

F 

T= 5 

Cette grandeur homogène à une force par unité de surface s’exprime en N m -2 (ou en Pa). 


8.1.3 Loi de Hooke 


Dans le domaine de l’élasticité linéaire, c’est-à-dire pour des tensions mécaniques et des 
déformations pas trop importantes, la déformation est proportionnelle à la contrainte. Cette 
propriété est exprimée par la loi de Hooke : 


ou de manière équivalente 


1 F 1 
£ “ Ë~S ~E T 


F „ 

T= s = E£ 


Dans cette expression la constante E, caractéristique du matériau constituant le fil, s’appelle 
le module d’Young; elle a les dimensions d’une force par unité de surface (en N m~ 2 ou en Pa). 


8.1.4 Coefficient de Poisson 

En plus de l’allongement selon Ox, il se produit un rétrécissement dans les directions per- 
pendiculaires à l’axe du fil. Si d est la dimension latérale du fil, on a : 

A d _d! -d _ AU 
~d ~ d ~~ V ~t 

Le coefficient v est appelé coefficient de Poisson ; c’est un nombre positif toujours inférieur à 
0.5. Pour les métaux usuels, sa valeur est v = 0.3. 


H. Djelouah 



8.2 Onde plane longitudinale 


69 


8.1.5 


Loi de Hooke pour les forces tangentielles 




Barreau soumis à une force de cisaillement. 

Si nous appliquons à la surface S d’un barreau une force F, non plus normale mais tangen- 
tielle, la longueur du barreau ne change pas, seules les arêtes normales au plan d’application de 
F tournent d’un angle a. 

Dans l’approximation de l’élasticité linéaire (cas des petites valeurs de a), l’angle a est 
proportionnel à la force appliquée. La loi de Hooke s’exprime dans ce cas par : 

^ = G a ~ G tan (a) = G ^ 

où G est le module de cisaillement ; il s’exprime également en N ■ m -2 . 

Le tableau ci-dessous donne les valeurs caractéristiques des modules de quelques matériaux 
courants. 


Matière 

p ( kg m 3 ) 

E ( N m" 2 ) 

G ( N m -2 ) 

Acier 

7.8xl0 3 

2.2xlO n 

0.9xl0 n 

Fer 

7.85 xlO 3 

2xlO n 

0.8xl0 n 

Aluminium 

2.7xl0 3 

0.7xl0 n 

0.26xl0 n 

Cuivre 

8.93xl0 3 

1.22xlO n 

0.42x 10 11 


Caractéristiques mécaniques de quelques métaux usuels. 


8.2 Onde plane longitudinale 

8.2.1 Equation de propagation 

Nous allons considérer ici uniquement le cas des ondes planes longitudinales se propageant 
dans un barreau. Considérons un élément du barreau de section S compris entre deux sections 
d’abscisses respectives x et x + Ax. La section d’abscisse x est soumise à une traction F ( x ) de 
la part de la partie gauche du barreau, tandis que la section d’abscisse x + Ax est soumise à une 
traction F (x + Ax) de la part de la partie droite du barreau. Sous l’action de ces deux forces 
les deux sections se déplacent respectivement de u x (x) et u x (x + Ax) le long de l’axe ( Ox ). Le 
déplacement de particules u x est une fonction à la fois de la coordonnée x et du temps t. Si p 
est la masse volumique du barreau, la masse de l’élément d’épaisseur Ax s’écrit : 

Am = pS Ax 
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Écrivons la relation fondamentale de la dynamique pour cet élément : 

F) 2 v 

p S Ax —Z = F (x + dx) - F (x) 
dt z 

Dans le cas où Ax est faible, la force F (x + Ax) agissant en x + Ax s’écrit : 


dF 

F (; x + Ax) = F (x) + Ax 

ax 

et la relation fondamentale de la dynamique devient 

ri d 2 u x dF 
pS Ax —y = — Ax 
dx z dx 

Sachant que la force normale F est reliée à la déformation e 

,du x 


F = SEe = SE- 


dx 


on obtient 


par la loi de Hooke : 


d 2 u x _ p_ d 2 u x _ 
dx 2 E dt 2 

On retrouve l’équation des ondes (équation de d’Alembert) qui décrit la propagation d’un 
ébranlement à la vitesse 



Les ondes décrites par cette équations correspondent à une grandeur physique vectorielle, 
le déplacement de particules u , orientée le long de la direction de propagation. Ce type de 
mouvement ondulatoire est dit longitudinal. 


8.2.2 Ondes progressives harmoniques 
Définition 

Dans le cas d’une onde progressive sinusoïdale, le déplacement de particules s’écrit en notation 
complexe : 


u x (x,t) = U 0 e i ^ t ~ k ^ 

où k = y = y es t l 0 module du vecteur d’onde, À étant la longueur d’onde. 

La composante F x de la force exercée en x par la partie gauche sur la partie droite est 

F x (x, t ) = -SE ^ = ikE u x (x, t) = ikEU 0 
dx 


Impédance 


On appelle impédance en un point le rapport de l’amplitude complexe de la force F x à 
l’amplitude complexe de la vitesse de particule ù x : 


Z (x) = 


Fx 

Ûx 


Dans le cas d’une onde progressive harmonique, on a : 

du x . . 

u x = u x (x,t) 
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d’où l’expression de l’impédance en un point 

Z(x) = SpV = Sy/pË 

On appelle impédance caractéristique du milieu constituant le barreau, la quantité 

Z c = yfpË = pV 

D’où l’expression de l’impédance mécanique du barreau : 

Z (x) = SZ C 

Nous remarquons que l’impédance caractéristique ne dépend pas de la position ; cette pro- 
priété est caractéristique des ondes planes progressives. Dans les autres cas, en présence de 
réflexion par exemple, l’impédance en un point dépend de la coordonnée x. Notons que l’impé- 
dance caractéristique du milieu Zc s’exprime en N m -2 tandis que l’impédance Z du barreau 
s’exprime en kg s -1 . 


8.2.3 Réflexion et transmission 

Réflexion et transmission entre deux barreaux semi-infinis 

Soit deux barreaux de même section S et de longueur semi-infinie, reliés en x = 0. Leurs 
masses volumiques sont respectivement p\ et p 2 ■ Lorsqu’une onde longitudinale venant de — oo 
se propage vers x = 0 dans le premier barreau , elle donne naissance au point de jonction, x = 0, 
à une onde réfléchie et une onde transmise. L’écriture de la continuité du déplacement et de la 
contrainte mécanique en x = 0 permet d’obtenir le coefficient de réflexion R u et le coefficient de 
transmission T u définis respectivement par : 


T u 


Ut 

Th 


où Ui , Ur et Ut sont les amplitudes des déplacements associés respectivement à l’onde 
incidente, l’onde réfléchie et l’onde transmise. On en déduit : 


R 


U 


Z\ — Z-2 
Z\ + Z 2 


T u 


2Zi 

Z\ + Z 2 


où Z\ et Z 2 sont les impédances respectives du premier et du deuxième barreau. 


Réflexion à l’extrémité d’un barreau terminé par une masse M. 

Considérons un barreau (p, S , E ) de longueur supposée semi-infinie (x < 0) terminé en son 
extrémité x = L par une masse M. Le barreau est le siège de la propagation d’un onde incidente 
et d’une onde réfléchie ; l’onde résultante s’écrit : 

u (x, t) = Ui + U R e^ t+kx) 


Écrivons la relation fondamentale de la dynamique pour la masse M : 


M 


d 2 u 

df 2 


= ES 


x=L 


du 

dx 


x=L 
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On en déduit l’expression du coefficient de réflexion au niveau de la masse M : 

Ur yfpE-itfr 

" U ■ æ=0 ^Æ + i^ 

On constate que le module du coefficient de réflexion \R U \ est égal à 1 ; il s’agit donc d’un 
phénomène de réflexion totale. On retrouve un résultat assez général qui est que le coefficient de 
réflexion est égal à la différence des impédances sur la somme des impédances et que dans le cas 
où l’impédance terminale est imaginaire, il se produit un phénomène de réflexion totale. 

8.2,4 Oscillations libres d’un barreau 

Considérons un barreau de longueur L dont l’une des extrémités (x = L) est libre et l’autre 
(x = 0) est fixée à un bâti rigide . Lorsque le barreau est soumis à une déformation longitudinale 
initiale, il est le siège d’ondes longitudinales se propageant dans le sens des x croissant et des x 
décroissants. Chacune de ces ondes subit une réflexion totale aux extrémités du barreau, donnant 
ainsi une multitude d’onde se propageant dans les deux sens. Ces différentes ondes interfèrent 
entre elles. 



Oscillations libres d’un barreau 


Pour ces raisons, nous recherchons une solution de l’équation d’onde sous la forme : 

u x (x, t) = A e A^~ kx ) + B e^ ut+kx) 

Dans cette dernière expression le vecteur d’onde k est donné par la relation de dispersion 
k = y tandis que A et B sont deux nombres complexes qui représentent l’amplitude complexe 
résultante de toutes les ondes se propageant respectivement dans le sens des x croissant et dans 
le sens des x décroissant. Ces deux nombres complexes dépendent des conditions aux limites en 
x = 0 et x = L , qui sont traduites par les équations suivantes : 

f u x (x = 0) = 0 

\ F(x = L) = SElfe | I=I = 0 

La première condition impose 

A + B = 0 


d’où B = —A, et 


u x (x, t) = A (e~ ikx - e +ikx ^j e iujt = -2iAsm (kx) e iuJt 
La seconde condition impose 


-2 ikSEA cos ( kL ) e iut = 0 


H. Djelouah 




8.2 Onde plane longitudinale 


73 


soit 

cos ( kL ) = 0 


dont les racines sont : 


< 2 " +1 >â 


n étant un nombre entier. Les pulsations correspondantes permettant de satisfaire les conditions 
aux frontières sont 


cc n — kn V — (2ti -f- 1 ) 


ttV 

~2L 


Ces pulsations "permises" sont appelées les pulsations propres. Il 3^ a donc une infinité de 
pulsations propres et la solution générale s’écrit : 


OO 

u x ( x , t ) = — 2iA n sin ( k n x ) e luJnt 

n = 0 

En revenant à la notation réelle, c’est-à-dire en prenant la partie réelle de l’expression précé- 
dente, on obtient : 


OO 

u x (*,i) = n cos (oj n ) + b n sin (uj n t)] sin ( k n x ) 

71=0 


où les coefficients a n et b n sont reliés aux nombres complexes A n par : —2 iA n = a n — ib n avec 

7 T ttV 

k n = (2n + 1) — et u n = k n V = (2 n + 1) — 

Les oj n sont les pulsations propres. Les coefficients a n et b n sont déterminés par les conditions 
initiales du mouvement. Supposons qu’à t = 0 nous imposions aux différents points du barreau 
un déplacement initial 

u(x, 0 ) = uq(x) 


et une vitesse initiale 


û (0, t) = v 0 ( t ) 

Dans ce cas nous aurons les conditions initiales suivantes : 


OO 

u o ( x ) = I] a n sin (k n x) 

n = 0 
00 

Vo ( X ) = Z u n b n sin (k n x) 

n = 0 


On doit inverser ces équations pour obtenir les coefficients a n et b n . La méthode de Fourier 
décrite dans le chapitre précédent, permet d’obtenir : 


2 f r 7txi 

a n = £ J Uo (x) sin |^(2n + 1) — j dx 

bn = [ V 0 (x) sin [(2n + 1) Z-] dx 




’o 


2L\ 
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8.2.5 Oscillations forcées d’un barreau de longueur finie 



Oscillations forcées d’un barreau de longueur finie. 


Considérons un barreau dont l’extrémité en x = 0 est fixée rigidement tandis que l’extrémité 
située en x = L est soumise à une force sinusoïdale F (t) = Fq cos (f U). 

En notation complexe les conditions aux frontières s’écrivent : 

r u x (o,t) = o 

1 F x (x = L,t) = F 0 e iQt 


Le milieu de propagation étant de longueur finie, le déplacement de particules u x ( x , t) s’écrit 
en régime permanent : 

u x ( x , t) = A e i{çlt ~ Kx) + B e i{ - m+Kx) 

avec K = y et 


du T 

F x = ES —F = -2 iKES 
ox 


g i(Çlt—Kx ) g i(flt+Kx ) 


En tenant compte des conditions aux limites on obtient les expressions de A et B : 

F 0 e iKL 


A = 
B = 


2 iKES cos (KL) 
F 0 e~ ŒL 

2iKES cos (KL) 


D’où l’expression du déplacement de particules : 


u ( x , t) = U ( x ) e lÇlt 


où l’amplitude U (x) s’exprime par : 

= Fo sin[ifæ] 

1 ; KES cos (KL) 

On obtient un phénomène d’ondes stationnaires dont les amplitudes et les positions des nœuds 
et des ventres pour le déplacement de particules sont donnés par le tableau suivant : 



Noeuds 

Ventres 

Position 

Xn = L - (2 n + 1) | 

x n = L — 

Amplitude 

t^min = 0 

TT — F 0 

u max — KES cos (KL) 


La résonance apparaît pour des fréquences d’excitations f î égales à l’une des pulsations 
propres du barreau Çî = oj n = (2n + l)^jr. Pour ces fréquences, l’amplitude t/ max au niveau 
des ventres d’élongation devient infinie. 
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8.3 Ondes élastiques transversales 



Nous allons analyser le problème des ondes élastiques transversales dans un barreau solide. 
Lorsque l’extrémité du barreau est soumise à une force de cisaillement, parallèle à la section S, 
nous pouvons supposer que chaque section du barreau se déplace de bas en haut et de haut en bas 
sans mouvement horizontal. Appelons u z le déplacement transversal d’une tranche d’épaisseur Ax 
à un instant donné. Le déplacement u z est fonction de la position x sinon il correspondrait à un 
déplacement de l’ensemble du barreau parallèlement à lui même. Il en résulte une déformation 
de cisaillement. Chaque tranche d’épaisseur Ax est soumise aux forces antagonistes F (x) et 
F (x + Ax) qui sont tangentes aux sections et qui sont produites par les portions du barreau qui 
sont situées de chaque côté du barreau. Il existe entre la force tangentielle de cisaillement et la 
déformation de cisaillement, une relation analogue à la loi de Hooke : 

F _ du z 
S du x 

où G est un coefficient caractéristique du matériau, que l’on nomme module de cisaillement. 

La force résultante sur la tranche Ax est : 

dF 

F (x + Ax) — F (x) = — Ax 

dx 

La relation fondamentale de la dynamique pour la tranche Ax s’écrit : 

fîF rfin 

F (x + Ax) — F (x) = — — Ax = pS Ax * 

ox dt z 

dF_ d 2 u z 

dx ^ dt 2 

La dérivation de la loi de Hooke, par rapport à x donne : 


— = ° 2uz 

dx dx 2 


Par substitution, on obtient 

d 2 u z _ p d 2 u z _ 
dx 2 G dt 2 

A nouveau nous obtenons l’équation de propagation d’une onde, de cisaillement cette fois, 
qui se propage à la vitesse 



La composante F z de la force exercée par la partie gauche du barreau sur l’élément d’épaisseur 
Ax est définie par : 


F z 


-SG 


du z 

dx 
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On peut montrer que l’impédance mécanique du barreau Z s’écrit : 

F ' 

Z = — = S z c 
u z 

où Z c = pV = \fpE est l’impédance caractéristique du matériau constituant le barreau. 

8.4 Modèle de la chaîne linéaire 

Dans ce paragraphe on développera sommairement le modèle de la chaîne linéaire d’atomes 
qui permet de décrire le phénomène de propagation des ondes élastiques dans les solides. Consi- 
dérons un barreau de section S et de longueur L "très grande", taillé dans un solide cristallin 
monoatomique. Si on applique sur la face située en x = 0 une force de compression F(t), une 
onde longitudinale va se propager le long de l’axe Ox. Étudions le phénomène qui se produit dans 
le barreau avant que le front d’onde n’ait atteint son extrémité. Nous supposons que le barreau 
est un solide monocristallin monoatomique et que les atomes le constituant sont disposés aux 
noeuds d’un réseau tridimensionnel régulier. En se propageant une onde plane longitudinale fait 
osciller simultanément tous les atomes se trouvant dans un plan perpendiculaire à la direction de 
propagation Ox. Puisqu’ils se déplacent tous en bloc, ces atomes sont équivalents à une masse 
unique M concentrée sur l’axe Ox. Chacun des plans transversaux d’atomes est maintenu en 
place par les forces de liaison qui agissent comme des ressorts placés en parallèle ; ils sont donc 
équivalents à un ressort unique de raideur K. Il suffit donc de raisonner sur une seule chaîne 
d’atomes identiques et équidistants, séparés à l’équilibre par une distance a. 

8.4.1 Modélisation microscopique du problème et mise en équations. 

Uq U\ U 2 U 3 1 U n ^n+l 


>o 


M 


x 0 


K 


Xi 


X2 


X3 


a 

>► 


^n— 1 




*£n+ 1 


Modèle de la chaîne linéaire 


X 


Le solide est constitué d’une chaîne comportant une infinité d’atomes, assimilés à des points 
matériels de même masse M , reliés par des ressorts identiques de longueur à vide a et de raideur 
K, et susceptibles de se déplacer sans frottement le long de l’axe Ox. 

Ces ressorts fictifs modélisent, dans l’approximation linéaire, les actions subies par les atomes 
lorsqu’ils se déplacent au voisinage de leurs position d’équilibre. C’est ce couplage entre les 
différents oscillateurs qui permet la propagation d’une onde élastique. On repère les positions à 
l’équilibre de l’atome d’ordre n par : x n \é qu = n a. Hors équilibre, les abscisses de ces atomes 
sont : x n = x\é qu + u n (t) = n a + u n (t) ; les déplacements particulaires u n (t) supposés faibles 
devant a, dépendent à la fois du temps et de la position à l’équilibre na. 

L’énergie potentielle de ce système s’écrit : 

U = ht (■ u 1 - U 0 ) 2 + { Ui - u 0 ) 2 H h ht (■ u n - u n - 1 ) 2 + (u n+1 - u n f H 

tandis que l’énergie cinétique s’écrit : 

T = -Mû § + - Mù\ -| h -Mùn H 
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D’où la fonction de Lagrange : 

C = T-U = Y j \mùI - \k K - U n - 1) 2 

n n 

L’écriture de l’équation de Lagrange permet d’obtenir l’équation différentielle du mouvement 
pour l’atome d’ordre n : 


M ^ — L (u n u n — 1 ) K (u n u n -\- 1 ) 

M'Ü n = -K (2 U n - Un- 1 - U n+ 1) 

8.4.2 Solution en régime permanent sinusoïdal 

Imposons à la particule d’ordre n = 0 une oscillation sinusoïdale : 

uq = ao e lut 


En régime permanent et en l’absence d’atténuation, le mouvement de chaque atome est 
identique à celui de la particule d’ordre n = 0 à un déphasage près qui correspond au temps de 
propagation de l’onde : 




i(iüt—kx n ) 


où k est le vecteur d’onde et x n = na. L’équation différentielle du mouvement s’écrit alors : 

(2 K - Mu 2 ) - K (e ika + e“ ifca ) = 0 


soit 


d’où 


(2 K - Mu 2 ) - 2/L cos (fca) = 0 


Mu z = AK sin z ( 


ou encore 



Cette dernière relation est appelée relation de dispersion. Sa représentation graphique est 
appelée diagramme de dispersion du milieu de propagation (figure ci-dessous). 
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Comme |sin (^f )| < 1, la fréquence des ondes élastiques pouvant se propager dans un cristal 
est limitée par la fréquence de coupure : 


f<fc 


Wç 
27 T 


i 

7T V M 


L’intervalle — ^ < k < ^ est appelé première zone de Brillouin. La longueur d’onde la plus 
courte que puisse transmettre le réseau est : 

, 2tt „ 

Xc — ~r — 2 a 


où k c = 

c a 

Nous remarquons que u n peut s’écrire également : 


iuj 

U n — CLQ ^ 



où v<f, = j représente la vitesse de phase. 

Pour les faibles valeurs du nombre d’onde ( ka « 1 ), la relation de dispersion s’écrit : 


c jü — 



C’est l’équation d’une droite dont la pente Vo = fly f représente la 
des ondes élastiques de basses fréquences. La fréquence de coupure peut 
Vo : 

A = ^ 

7 ra 


vitesse de propagation 
s’écrire en fonction de 


8.4.3 L’approximation d’un milieu continu. 

Dans un réseau cristallin, la distance a entre deux atomes est typiquement de l’ordre de 
lCT 10 m, distance très inférieure aux dimensions caractéristiques des phénomènes de propagation 
à étudier (qui sont plutôt de l’ordre de grandeur du mm à quelques mètres). Ceci suggère de 
considérer la chaîne d’atomes comme un milieu continu, dont les déplacements sont décrits par 
une fonction continue u(x,t) prenant les valeurs u n {t) aux points d’abscisses x n = na. 

L’équation du mouvement de la particule d’ordre n peut se réécrire sous la forme : 

MÜ n - I< (u n+ 1 - Un) + I< ( U n - U n+ i) = 0 

Si À >> a, la différence u n + 1 — u n est un infiniment petit qui, au premier ordre près, peut 
s’écrire : 


^n+1 ^n+1 ^ n 

^n+1 %n & 

U n U n —1 Un "Un — 1 

•En %n—l U 



# 71+1 



#71 


Dans ces conditions, l’équation différentielle du mouvement devient : 


i. T d 2 u 

M w = Ka 




n 
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Au premier ordre près, on peut considérer que : 


\OX J x n +l \OXJxn 
^n+1 %n 


d 2 u 
dx 2 


et l’équation du mouvement s’écrit alors : 

d 2 u Ka 2 d 2 u 

dt, 2 M dx 2 

Cette équation est équivalente à l’équation de propagation : 

d 2 u d 2 u 
dt 2 0 dx 2 

où Vq = est la vitesse de propagation des ondes élastiques de basse fréquence. 

Ordres de grandeur : La vitesse Vo des ondes élastiques de basse fréquence dans les solides est 
en général comprise entre 1000 m s -1 et 10000 m s -1 . Les distances entre atomes sont de quelques 
angstroms. Avec Vo = 5000 m s -1 et a = 5À, la fréquence de coupure est f c = 3.2 x 10 12 Hz. 
Cette valeur est si élevée que pour le domaine des fréquences inférieures à quelques gigahertz, on 
se situe au tout début de la courbe de dispersion, là où la pulsation est proportionnelle au nombre 
d’onde. La longueur d’onde À = v jt , comprise dans ces cas entre quelques mm et quelques gm, 
est très grande devant les distances interatomiques, si bien que pour l’onde le milieu apparaît 
comme continu. 


8.5 Exercices 


Exercice 1 : 


X x+Ax 

[ j I ^ section S 


L ~>‘ 

U(x,t) 


Soit un barreau de grande longueur et de faible section S, soumis à des vibrations longitudinales. 
La masse volumique du barreau et son module d’Young sont respectivement p et E. Considérons 
un petit élément d’épaisseur faible égale à Ax, situé à la position x. Lorsque cet élément est 
atteint par la vibration, il se met en mouvement sous l’action de la résultante des forces exercées 
par les parties du barreau qui se trouvent à gauche et à droite,. En raison de l’élasticité du 
milieu, les deux sections situées respectivement en x et x + Ax subissent des déplacements 
différents notés respectivement u(x, t ) et u(x + Ax, t) ; l’élément subit donc une déformation en 
plus du mouvement de translation autour de sa position d’équilibre. On négligera dans la suite 
la déformation latérale. 

1. Montrer que u(x,t ) correspond à une onde longitudinale qui se propage à une vitesse V 
que l’on précisera. Dans le cas d’une onde progressive sinusoïdale, calculer l’impédance en 
un point quelconque du barreau. 

2. Ce barreau, dont la longueur est égale à L, est encastré dans une paroi rigide. Son autre 
extrémité est libre. 
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F(t) ) 




X 

v 


0 

L 


(a) Oscillations libres : Ecrire l’expression générale de u(x,t) sous la forme d’ondes si- 
nusoïdales de pulsation u. Ecrire les conditions aux limites en x = 0 et x = L. En 
déduire les pulsations propres du barreau. 

(b) Oscillations forcées : L’extrémité du barreau en x = 0 est soumise à une force hori- 
zontale F(t) = F(j cos (fit) . Montrer que le déplacement en chaque point s’écrit sous 
la forme u(x,t ) = U(x)cos(Qt + <j>) ; quelle est l’expression de l’amplitude U(x) et de 
la phase 4> . Quelle est la position des minima et des maxima de vibration ? Quelles 
sont les amplitudes U min et Umax en ces points. Pour quelles pulsations obtient-on le 
phénomène de résonance ? 

Exercice 2 : Un barreau de longueur L est fixé rigidement en x = 0 ; son extrémité en x = L 
est libre. 

1. Montrer que seuls les harmoniques impairs sont permis. 

2. Déterminer la fréquence fondamentale du barreau s’il est constitué d’acier et que sa lon- 
gueur est L = 0.5 m. 

3. Si une force statique est appliquée à l’extrémité libre du barreau de telle sorte à déplacer 
cette extrémité d’une longueur égale à h, montrer que, lorsque le barreau vibre suite à une 
suppression brusque de cette force, les amplitudes des différents harmoniques de vibration 
sont donnés par A n = [8h/nn 2 ] sin(n7r/2) . 

4. Déterminer ces amplitudes si la force statique vaut Fs = 5000 N et si la section droite du 
barreau a une surface égale à s = 0.00005 m 2 . 

N. B : Rechercher préalablement les caractéristiques mécaniques de l’acier dans la bibliographie. 

Exercice 3 : Un barreau d’acier de 0.0001 m 2 de section droite et de 0.25 m de longueur est 
libre de se mouvoir en x = 0 et il est chargé à l’autre extrémité , en x = 0.25 m, par une masse 
de 0.15 kg. 

1. Calculer la fréquence fondamentale des vibrations longitudinales de ce barreau. 

2. Déterminer la position à laquelle ce barreau doit être fixé pour causer la plus faible per- 
turbation de son mode fondamental de vibration. 

3. Quand ce barreau vibre dans son mode fondamental, quel est le rapport de l’amplitude du 
déplacement de l’extrémité libre à celle qui est chargée par la masse? 

4. Quelle est la fréquence du premier harmonique de ce barreau ? 

Exercice 4 : Une masse m = 2 kg est suspendue à un fil d’acier de longueur £ = 1.0 m et de 
section droite s = 0.00001 m 2 . 

1. Calculer la fréquence fondamentale des oscillations verticales de la masse en la considérant 
comme un simple oscillateur. 

2. Calculer la fréquence fondamentale des oscillations verticales du système constitué par un 
barreau vibrant fixé à une extrémité et chargé par une masse à l’autre extrémité. 

3. Montrer que pour kL < 0.2, l’équation obtenue à la question précédente peut se mettre 
sous la forme : u;q = W s/m. 
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Exercice 5 : Soit une chaîne linéaire de particule identiques de masse rn. A l’équilibre, la masse 
d’ordre n est située, sur l’axe Ox, au point d’abscisse x n = n ■ a, x n étant négatif ou nul (Figure 
a). Les masses sont reliées entre elles par des ressorts identiques de raideur K, dont la longueur 
au repos est a. Soit u n l’écart de la particule d’ordre n par rapport à la position d’équilibre. 




1. Ecrire l’équation du mouvement de la particule d’ordre n. 

2. En cherchant les solutions sous la forme d’ondes progressives sinusoïdales planes telles que 

u n = où jj q es t un e constante réelle, déterminer la relation de dispersion 

UJ = tjü(k). 

3. La particule en xo = 0 est fixe. On envoie une onde incidente de la gauche vers la droite 
décrite par u n = u^^t.-kna) _ Déterminer l’amplitude U' et la phase 4> de l’onde réfléchie. 

4. La particule en x = 0 est maintenant liée à un bâti fixe par l’intermédiaire d’un ressort 
de raideur K' et d’un amortisseur de coefficient de frottement visqueux a (Figure b). On 
choisit K' et a pour qu’il n’y ait pas d’onde réfléchie. 

(a) Ecrire l’équation du mouvement de la masse située en x = 0. 

(b) On envoie une onde incidente comme précédemment. Déterminer les coefficients I\' 
et a pour que cette masse vibre sous l’effet de la seule onde incidente (pas d’onde 
réfléchie). 


Exercice 6 : Soit une corde plombée telle que représentée sur la figure ci-dessous, où x n 

représente la position de la masse m et u n son déplacement. Les plombs sont espacés l’un de 
l’autre de a et reliés entre eux par des segments de fil de tension Tq que l’on supposera constante 
lors des déplacements des masses. Cette corde plombée de longueur L telle que L » a, subit 
une vibration transversale. 



Dans le cas des vibrations de faible amplitude, on demande de : 

1. Ecrire l’équation du mouvement du plomb d’ordre n. 

2. Chercher les solutions de cette équation sous la forme d’une onde progressive. 

3. En déduire la relation entre la pulsation oj et le vecteur d’onde k. 

4. Montrer que l’on peut limiter l’étude à l’intervalle 0 < \k\ < ir/a. Tracer la courbe t o = oj(k). 

5. Calculer la pulsation de coupure lom ■ 
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6. Définir la vitesse de phase et la vitesse de groupe Vq- Décrire le mouvement lorsque 
k = ir/a. 

7. Dans le cas des très grandes longueurs d’onde, que deviennent V# et Vq ? 
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Chapitre 9 


Ondes acoustiques dans les fluides 


9.1 Introduction 

Les ondes acoustiques sont des ondes élastiques qui se propagent dans les fluides (gaz ou 
liquides). Il est donc possible d’obtenir l’équation d’onde qui régit la propagation des ondes 
planes dans un fluide par la même démarche que celle que nous avons utilisée pour établir 
l’équation de propagation des ondes transversales dans une corde. 

Dans la suite, nous utiliserons les symboles suivants pour étudier l’onde acoustique qui se 
propage suivant l’axe des x : 

x : coordonnée à l’équilibre d’une particule du milieu. 

u x : composante suivant l’axe des x du déplacement de particule par rapport à la position 
d’équilibre. 

Pc : masse volumique du fluide à l’équilibre 

P : pression instantanée en un point quelconque 

Po : pression à l’équilibre 

p = P — Po : surpression ou pression acoustique 

c : vitesse de propagation de l’onde 

On entend par particule, un élément de volume contenant des millions de molécules de telle 
sorte qu’il puisse être considéré comme continu, mais toutefois suffisamment petit pour que les 
grandeurs acoustiques comme la pression, la masse volumique et la vitesse de particule puissent 
être considérées comme constantes dans cet élément de volume. Dans ce qui suit, nous négligerons 
les effets de la gravitation de telle sorte que Po et po sont uniformes dans tout le milieu. On 
suppose d’autre part que le milieu est homogène, isotrope et parfaitement élastique, c’est-à-dire 
non dissipatif. 


9.2 Equation d’onde 


Considérons le cas d’une onde plane émise dans un fluide par une membrane vibrante plane. 
Lorsque celle-ci est au repos, la pression dans le fluide est uniforme et égale à Pq. En se déplaçant, 
par exemple dans le sens des x positifs, la membrane comprime la couche de fluide adjacente. 
Cette situation est instable : le fluide se détend en comprimant à son tour la tranche voisine. 
L’onde progresse ainsi de proche en proche par une succession de compressions et de détentes. 
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P(x) . 

1 

1 

1 

1 

1 

1 

1 

1 

1 

.P(x+dx) ! 

f] 

S 


1 

1 

X 


x+dx 



u x (x) 

u x (x+dx) 



dx 

Propagation d’une onde acoustique 


Soit une tranche de fluide de petite épaisseur Ax située à l’abscisse x. Lorsque la perturbation 
atteint ce point, les forces agissant sur cette tranche ne s’équilibrent plus et elle se met en 
mouvement. Soit u x (x,t ) le déplacement à l’instant t du plan d’abscisse x. Soit F x (x,t ) et 
F x (x + Ax, t ) les forces agissant sur la tranche de fluide respectivement en x et x + Ax. Ces 
forces s’expriment par : 


F x ( x , t) = S P (x, t ) 

F x (x + Ax, t) = —S P (x + Ax, t ) 

La résultante de ces deux forces est : 


A F x = F x ( x , t) + F x (x + Ax, t ) 

AF X = —S [P (x + Ax, t) — P (x, t)] 

En faisant un développement en série de Taylor au premier ordre de P (x,t), on obtient : 

dP 

P(x + Ax, t) = P (x, t ) + — — Ax 

ax 

D’où : 


dp 

A F x = -S — Ax 
ax 

Comme P = Pq + p, la force résultante s’exprime par : 

A F x = -S ^ Ax 
ax 

Sous l’action de cette force, la tranche de fluide subit une accélération et en écrivant la relation 
fondamentale de la dynamique, on obtient : 

d 2 u x dp 

Am — = -S — Ax 
at. z ax 

Le fluide étant compressible, le déplacement du plan d’abscisse x + Ax est différent du 
déplacement du plan d’abscisse x et il vaut u x (x + Ax,t). De nouveau un développement en 
série de Taylor au premier ordre permet d’écrire : 

dux 

u x ( x + Ax, t) = u x (x, t ) H — — Ax 

ax 

Pour prendre en compte la compressibilité du fluide, calculons la dilatation volumique subie 
par la tranche de fluide . Soit Av$ = S Ax, le volume à l’équilibre et soit An, le volume en cours 
de mouvement, avec : 
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Ax 

Av 

Av 

Av 


S [x + Ax + u x (x + Ax) — x — u x (x) 
S [Ax + u x (x + Ax) — u x (x)] 


5 


* du x . 
Ax + — — Ax 
ax 


Ax 0 + < ^Ax 0 
ax 


On en déduit la dilatation volumique 


Av — Axo 
Ax 0 


dx 

Rappelons que pour un fluide compressible, la surpression p est reliée à la dilatation volumique 
6 par la relation 

p = —k 9 

où k est le module de compressibilité. On obtient ainsi : 


du x 


Remarque : On utilise souvent le coefficient de compressibilité x 
Les deux équations 


1 

K 


PO 


d 2 u x 

dt 2 


dp 

dx 

du r 


p = —K 


dx 


constituent les deux équations fondamentales de l’acoustique. 

En remplaçant dans la première équation p par son expression tirée de la seconde équation 
on obtient l’équation de propagation : 


d 2 p 1 d 2 p _ 
dx 2 V 2 dt 2 

où V représente la vitesse de propagation définie par 


V = 



1 

ffipoX 


9.3 Vitesse du son 

Le phénomène de propagation étant un processus adiabatique, la relation liant la pression et 
le volume est 

Pv 7 = constante 

En calculant la différentielle, on obtient : 

x 7 dP + 7 -P x 7_1 dv = 0 
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Si l’on considère que dP représente la variation de pression au voisinage de la pression à 
l’équilibre Pq, on obtient : 

VqP + 7P0 1 An = 0 


d’où : 


An 

P = -7-n) — 


En tenant compte de la définition du module de compressibilité, on obtient : 

K=-= 7-P 0 
X 


D’où la vitesse du son dans un fluide : 


V = 



Exemple : Dans l’air, dans les conditions normales T = 20 °C et Po = 10 5 iV • m~ 2 , 7 = 1.4 et 
po = 1.29 kg ■ m~ 3 , on en déduit V — 330 m ■ s -1 

La valeur de la pression à l’équilibre dépend fortement de la température. Pour une mole de 
gaz parfait, on a : 


P 0 v 0 = RT 


d’où 


et 


Po = 


RT 

vo 



Le produit pqvq représente la masse molaire M du gaz ; d’où : 


V = 



Dans un gaz parfait, la vitesse de propagation du son est proportionnelle à la racine carrée 
de la température mesurée en kelvin (K). 


9.4 Onde progressive sinusoïdale 

9.4.1 Définition 

Une onde acoustique sinusoïdale, s’écrit : 


p(x,t ) = po cos | U) [ t — 
On définit le module du vecteur d’onde k par 


h — — 

k ~v 


x 


V)\ 


d’où 


p ( x , t) = po cos (ut — kx ) 
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En notation complexe, l’onde progressive sinusoïdale s’écrit 

p(x,t)=p 0 e i ^ t - kx ï 

La relation liant la pression acoustique et la compressibilité, à savoir 

du 

p = -k — 

ox 

permet d’écrire 

u(x,t ) = p(x,t)dx 

K J 

u{x,t) = ~J Poe i{wt ~ kx) dx 

u fat) = -y 2 - e i{üjt ~ kx) 

IKK 

u fat) = 

La dérivation de cette dernière expression par rapport au temps permet d’obtenir la vitesse 
de particules : 

*<*•*> = £ = 

On constate que pour une onde progressive la vitesse de particules est en phase avec la 
pression acoustique. 


9.4.2 Impédance acoustique 

On appelle impédance acoustique en un point le rapport de l’amplitude complexe de la 
pression à l’ampliude complexe de la vitesse de particule 

z( *) = ? 

Dans le cas d’une onde progressive, on obtient : 

Z(x) = Po V 

Le produit poV définit l’impédance acoustique caractéristique du fluide 

Z c = PoV 

On obtient une propriété de l’onde plane progressive : 

z fa) = z c y x 


9.4.3 Energie acoustique 
Densité d’énergie cinétique 

Soit un petit élément de volume vç, dont le dépélacement est u (x, t ) et dont la vitesse est 
ù ( x , t) ; il possède une énergie cinétique 

p 1 -2 
E c = ~Pov 0 u 


On définit l’énergie cinétique par unité de volume ou densité d’énergie cinétique 

„ E c 1 2 

£c = — = X POU 2 

vq 2 
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Densité d’énergie potentielle 

Soit un petit élément de volume vq. Sous l’action de la surpression p, cet élément se comprime 
ou se dilate en raison de la compressibilité du fluide. L’énergie potentielle emmagasinée est égale 
au travail fourni par la pression pour comprimer ou dilater le volume vq : 


Sachant que 


on en déduit que : 



v - v 0 
p = —k 


dv = Vq dp 

K 

d’où : 


E p = —[ p dp 

k J o 


TP V ° 2 

EP= T K P 

On en déduit la densité d’énergie potentielle : 


g = Ep = _t_ 2 

P Vq 2k 


Densité d’énergie 

La densité d’énergie est égale à la somme de la densité d’énergie cinétique et de la densité 
d’énergie potentielle 


£ — £ c -\- £ n 


„ 1 o 1 2 

£= 2' , »“ + ï?' 

Dans le cas particulier d’une onde plane progressive sinusoïdale, 

1 


£c ~ £ p- 2po Vi 


po cos 2 (eut — kx) 


et 


£ = — —k po cos 2 (ut — kx) 

PoV - 


On définit la valeur moyenne temporelle de la densité d’énergie comme 

1 [ t+T 

( £ ) = f J t £dt 

m 2 t r 

où T = — . 

UJ 
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On définit également la moyenne spatiale de la densité d’énergie : 

rx-\- À 


1 r 


8 dx 


Dans le cas d’une onde progressive, ces deux valeurs sont égales et valent : 


<£> = 


Po 


2 Po V 2 


Intensité 

On appelle intensité de l’onde acoustique la puissance qui traverse, par unité de temps, une 
surface unité perpendiculaire à la direction de propagation. 

Pour calculer l’intensité de l’onde calculons l’énergie qui traverse pendant un intervalle de 
temps une surface S perpendiculaire à la direction de propagation. 



► 

Vdt 

Flux de puissance 


Cette énergie dE est égale à l’énergie contenue dans un volume SV dt et elle est égale à 

dE = S S V dt 


D’où la puissance V traversant cette surface 

V = ^ =SSV 
dt 

On en déduit l’expression de l’intensité de l’onde acoustique 

lit) =8 V 

I (t) = cos 2 (ut — kx ) 

PqV 

On appelle intensité de l’onde acoustique la valeur moyenne 


'-è 


rt+T 


I ( t ) dt 


I = 


I = 


pI 

2 PqV 

A 

2Z r 
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Niveau sonore 

On définit le niveau sonore en décibels par 

N dB = 10 log ( I - 

Io est une intensité de référence correspondant à Iq = 10~ 12 W ■ m~ 2 . 

Exemple : Pour une fréquence / = 1kHz, le seuil d’audition est égal à 0 dB et le seuil de 
douleur est égal à 130 dB. Pour calculer l’intensité, l’amplitude de pression po la vitesse de 
particules et le déplacement de particules, dans le cas où V = 330 m/s et Z c = 411 rayleighs, on 
peut utiliser les relations suivantes : 

I = I o 10 O1N dB 

Po = \J2Z C I 
Po 


u = 


Z t 


u = 


2vr/ 


Pour chacun de ces deux cas, on obtient : 


NdB 

I ( W/m 2 ) 

Po ( Pa ) 

û ( m/s ) 

u ( m ) 

0 dB 

io - 12 

2.9 x IO" 5 

7 X IO" 8 

1. 1 x 10“ n 

130 dB 

10.0 

91 

0.22 

3.5 x IO" 5 


9.5 Réflexion- Transmission 

Onde incidente Onde transmise 

i =^> 

~ c5 X 

Onde réfléchie 

Pl’ V t P2.V 2 

Réflexion à un interface fluide-fluide 

Soit deux milieux fluides semi-infinis séparés par un surface plane. Choisissons un repère 
orthonormé de telle sorte que le plan yOz coïncide avec la surface de séparation. Lorsque une 
onde acoustique provenant de — oo, se propageant dans le premier dans la direction de l’axe des 
x arrive à la surface de séparation, elle donne naissance à deux ondes 

- une onde réfléchie qui se propage dans le premier milieu dans le sens des x décroissants. 

- une onde transmise qui se propage dans le second milieu dans le sens des x croissants. 

L’onde résultante dans le premier milieu (x < 0) est caractérisée par : 


pi (x, t.) = pi ( x , t) + p R ( x , t) 
pi (x,t) = P i e i ^ t ~ k ^ + P R e 


ûi (x,t) 


1 


p i{u>t—kix) 


i(uit+k\x) 

Pr 


Dans le deuxième milieu, on a 
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P2 {X, t.) = p T (• X , t ) 

P 2 (x,t) = P r e i ( ut - haX ï 

ù 2 (x,t) = _Lp Te *M-^) 
£-2 

Les relations de continuité à l’interface s’écrivent 


Pi (0 ,t) =p 2 (0 ,t) 
ûi (0, t) = Ù2 (0, t) 


On en déduit 


ou encore 


Pi + P R = Pt 
Y l (■ Pi ~ Pr ) = ~^ Pt 


Pr Pt 

1 H — — = — - 

Pi p u 

^ _ Pr _ Zi 

Pi ~ Z2^I 


On définit 

- le coefficient de réflexion pour la pression 




le coefficient de transmission pour la pression 


T P = 


Les deux relations de continuité s’écrivent alors 


Pt 

Pi 


1 + Rp = Tp 
1 -Rp = %T P 


On en déduit les coefficients de réflexion et de transmission 


Rp 


T P 


Z 2 — Z\ 
Z 2 + Z\ 
2Z 2 

Z 2 + Z\ 


En tenant compte des relations = Z\ ù*. pp = —Z\ ùp et pp = Z 2 ûp, on peut calculer les 
coefficients de réflexion et de transmission pour la vitesse de particules et pour le déplacement 
de particules : 


R 


u 


■u 


Ru 


Tu 


Z\ — Z 2 
Z\ + z 2 

‘ÏZ\ 

Z\ + z 2 
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En tenant compte des relations = P? /2Zi, I R = P^/2Zi et I R = P^/2Z-2, on peut calcul el- 
les coefficients de réflexion et de transmission pour l’intensité acoustique : 



It 



Z\ — Z 2 
Z\ + Z^_ 
4 Z\Z-2 
[Zi + z 2 ] 2 


9.6 Exercices 

Exercice 1 : Dans les conditions normales de température et de pression, l’air est caractérisé 
par une masse volumique à l’équilibre p = 1.21 kg/m 3 et une valeur de 7 = c p /c v = 1.402. 
Calculer la vitesse de propagation du son dans l’air et son impédance acoustique spécifique dans 
ces conditions de température et de pression (T = 20 °C et Po = 10 5 N/m 2 ). 

Exercice 2 : Soit une onde acoustique plane de fréquence / = 24 kHz, qui se propage dans 
l’eau avec une vitesse c = 1480 m/s. Elle véhicule une puissance moyenne de V = 1 W unifor- 
mément répartie sur une section circulaire de de diamètre D = 40 cm, normale à la direction de 
propagation. La fréquence de l’onde est / = 24 kHz. 

1. Calculer l’intensité acoustique; quel est en dB le niveau de l’intensité acoustique relative- 
ment à un niveau de référence Iq = 10 ~ 12 W/m 2 ? 

2. Calculer l’amplitude de la pression acoustique, l’amplitude de la vitesse de particules et 
l’amplitude du déplacement de particules. 

3. Comparer aux résultats que l’on aurait obtenus si cette onde se propageait dans l’air. 

Exercice 3 : Une onde acoustique plane se propageant dans l’eau arrive en incidence normale 
à la surface de séparation avec l’air. Calculer les valeurs numériques des rapports suivants : 
PrJ Pi, P T / PiiU R /UiiU T /Ui,Ü R /Ùi,Ü T /Ùi, I R / IiJ T / Ii où U, Ù, P et / représentent respecti- 
vement le déplacement des particules, la vitesse des particules, la pression acoustique et l’in- 
tensité acoustique. Les indices i, R et T se rapportent respectivement à l’onde incidente, l’onde 
réfléchie et l’onde transmise. 

Exercice 4 : Répondre aux mêmes questions que pour l’exercice précédent en supposant que 
Fonde acoustique se propage initialement dans l’air et se transmet dans l’eau. Comparer aux 
résultats de l’exercice précédent. 

Exercice 5 : Un tuj^au cylindrique de section S constante est rempli d’un gaz, de masse volu- 
mique p , où les ondes acoustiques peuvent se propager à la vitesse V. A l’une des extrémités 
du tuyau, en x = 0, on place un piston plan qui est animé, suivant l’axe Ox du tuyau, d’un 
déplacement sinusoïdal d’amplitude JJq et de pulsation ui . A la position x = L, le tuyau est 
fermé par une paroi rigide. 

1. Montrer que le déplacement peut se mettre sous la forme : u(x,t) = U(x) e l ^ ut+ ^ où U(x) 
est réel. Donner l’expression de l’amplitude U{x) et de la phase </> de Fonde résultante. 

2. En déduire les positions et les valeurs des maxima et des minima de l’amplitude U(x) en 
valeur absolue. 

3. Pour quelles fréquences obtient-on un phénomène de résonance ? Quelle est alors l’ampli- 
tude du déplacement au niveau des ventres de déplacement. 

4. Dans le cas où L = 13À/12, A étant la longueur d’onde, tracer sur un graphe les variations 
de \U{x)\ en fonction de x. Echelle : 1 cm = À/12. Quel est le nombre de maxima (ventres) 
et de minima (noeuds) ? 
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Exercice 6 : Une onde acoustique plane, sinusoïdale de pulsation eu , d’amplitude A\ . se propage 
suivant la direction Ox dans un milieu fluide de masse volumique p \ . Elle arrive sous incidence 
normale sur un second milieu fluide de masse volumique p 2 et d’épaisseur L 2 , déposé sur un solide 
parfaitement rigide (voir figure). On notera V\ et V 2 les vitesses de propagation respectives dans 
chacun des deux milieux fluides. Les vecteurs d’onde correspondants seront notés respectivement 
k 1 et k 2 - 


X' 

Milieu 1 

0 

Onde 

Incidente 

1 r 

Milieu 2 

i 

y 

L-2 

F 

Solide rigide 



f 


Dans chacun des deux milieux, la pression acoustique s’écrit respectivement 

pi (x, t ) = Ax + A[ e i ( wt + k i x ) 

P2 (; x , t ) = A 2 + A' 2 e^ t+k2x) 

1. Ecrire les expressions respectives, û\ (x,t) et Ù 2 (x,t) , des vitesses de particule dans le 
milieu 1 et dans le milieu 2. 

2. Ecrire les relations de continuité pour la pression acoustique et pour la vitesse de particule 
en x = 0. En déduire deux relations entre les coefficients A\ . A \ , A 2 et A' 2 . 

3. Ecrire la relation de continuité pour la vitesse de particule en x = L 2 . En déduire une 
relation entre les coefficients A 2 et A! 2 

A! 

4. En déduire le coefficient de réflexion R = En donner le module et la phase. 

5. Mesure de la vitesse V 2 

(a) Quelle est la condition pour que, dans le milieu 2, la pression possède un nœud en 
x = 0 et un seul ventre ? 

(b) Sachant que le premier nœud de pression dans le milieu 1 se situe en x = —L\, calculer 
la vitesse de propagation V 2 en fonction de V\, L\ et L 2 . 

(c) L’expérience est réalisée, dans les conditions ci-dessus, avec de l’eau comme milieu 1 
( V\ = 1500 m/s) et de la glycérine comme milieu 2. On mesure L 1 = 7.5 mm et L 2 = 
4.95 mm. Quelle est la vitesse de propagation V 2 des ondes acoustiques dans la glycé- 
rine ? 
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Annexe A 


Equations différentielles 


A.l Introduction 

Les oscillations linéaires des systèmes à un degré de liberté sont régies par des équations 
différentielles du second ordre à coefficients constants. Une équation différentielle linéaire du 
second ordre à coefficients constants est une relation entre une variable dépendante y et ses 
dérivées première et seconde par rapport à une variable indépendante t. qui peut s’écrire sous la 
forme : 


d 2 y , nx dy 2 _ ... 

^ + 2 Ô dt +U;oy - A{t) 

Les coefficients ô et oüq sont des constantes réelles positives. Dans les problèmes de vibration 
le paramètre t représente le temps et on note par convention : 


dy 

-dt =V 
d 2 y .. 

-dï =y 

D’où l’écriture de l’équation différentielle : 

ÿ + 2ôÿ + uly = A(t) 

A. 2 Equation homogène 

L’équation différentielle est dite homogène si A(t) = 0 : 

ÿ + 25ÿ + uoy = 0 

On recherche des solutions qui sont de la forme y(t) = e st . Dans ce cas l’équation différentielle 
s’écrit : 

[s 2 + 25 s + Wq] e st = 0 

Cette équation doit être vraie quel que soit t, ce qui implique : 

s 2 + 25 s + Wq =0 

On obtient ainsi une équation du second degré en s dite équation caractéristique. Les racines 
de cette équation caractéristique sont : 
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si = —5 + ^ 5 2 - ojI 



La solution de l’équation différentielle s’écrit alors : 


y(t) = A ie sit + A 2 e S2t 

A\ et A 2 sont deux constantes d’intégration que l’on peut déterminer à partir des conditions 
initiales : 


y (t = 0) = y 0 
ÿ(t = 0) = ÿ 0 

Sachant que ô et wo sont des nombres réels positifs, si et s 2 sont négatives ou complexes avec 
une partie réelle négative. La nature des solutions si et s 2 de l’équation caractéristique dépend 
de la valeur relative de <5 par rapport à loq . Ainsi trois cas sont à envisager : 

A. 2.1 Régime fortement amorti ( ô > cu 0 ) 

Dans ce cas, les deux racines si et s 2 sont réelles et négatives. L’écriture des conditions 
initiales donne un système de deux équations : 

Ai + A 2 = y 0 
siAi + s 2 A 2 = ÿ 0 

dont les solutions sont les constantes d’intégration A i et A 2 : 

, s 2 y 0 - ÿ 0 

Ai — 

«2 - Si 

A 2 = ÿ°~ sm 

S2 - Si 

y(t) est une fonction décroissant sans oscillations vers zéro lorsque t augmente. La forme 
exacte de y(t) dépend des valeurs de Ai et A 2 qui sont déterminées par les conditions initiales. 
Pour les conditions initiales particulières suivantes :yo / 0; ÿo = 0, la figure suivante représente 
les variations de y au cours du temps 
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A. 2. 2 Régime critique (S — ujq ) 

L’équation caractéristique possède une racine double : 

S\= S2 = -S 

La solution générale de l’équation différentielle est de la forme : 

y(t) = (Ai + A 2 t) e~ 5t 

Les constantes d’intégration A\ et A 2 sont déterminées à partir des conditions initiales et 
valent : 


Ai = yo 
A 2 = ÿo + S y 0 

y(t) et une fonction décroissant vers zéro quand t augmente. Il est aisé de vérifier que cette 
situation correspond à la décroissance la plus rapide de y(t). Ce cas correspond au régime critique. 
Pour le cas particulier de conditions initiales :yo / O; ÿo = 0, la solution est : 

y(t) = yoe~ 5t (l + ôt) 

La figure ci-dessous représente les variations de y(t). 



A. 2. 3 Régime pseudo-périodique ( ô < co 0 ) 

Dans ce cas si et S 2 sont complexes conjugués : 


si = — <5 + i\/ujQ - ô 2 


s 2 = —S — - ô 2 

où i représente le nombre imaginaire pur vérifiant la relation i 2 = —1. Posons 


ujd = V^o 


- ./,, 2_ ô 2 


La solution générale s’écrit alors : 
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y(t) = e~ st [Ai e iUDt + A 2 e~ iu}Dt ] 

Sachant que : 


e uo D t _ cos (w D t) + i sin (ojot) 
e -iu D t _ cos (uj D t) — i sin(ujDt) 

y(t ) s’écrit : 

y(t) = e~ St [(Ai + A 2 ) cos(o Jot) + i (Ai - A 2 ) sin(a;_Dt)] 

■y(t) étant réelle, les nombres complexes A\ et A 2 doivent vérifier les relations : 

Ai + A 2 : réel 
Ai — A 2 : imaginaire 

Ai et A 2 sont donc complexes conjugués et peuvent se mettre sous la forme : 


Ai = A' e i<fi 


A 2 = A' e~ i(p 


y(t ) s’écrit alors : 


y(t) 


A' 


-st 


e i(u D t+ip) _|_ e -i{tü D t+ip) 


soit finalement, en posant A = 2 A 1 : 


y(t) = A e St cos(o i£>t + (p) 

y(t ) peut être interprétée comme une fonction périodique de pulsation uip , de phase initiale 
(p et d’amplitude Ae~ 5t décroissant exponentiellement au cours du temps. On peut définir une 
pseudo-période 



lüd 


A et p) sont deux constantes d’intégration déünies à partir des conditions initiales : 


yo = A cos (ip) 

ÿo = -A cos (ip) [5 + c od tan(</?)] 

D’où l’on tire : 
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p = — arctan 


yo + àyo 

UDV 0 


4 = \/ {up yo) 2 + (j /0 + Syoÿ 

iüD 

La figure ci-dessous représente les variations de y(t) dans le cas particulier de conditions 
initiales :î/q 7 ^ 0; ÿo = 0. 



Remarque : Cas particulier où 5 = 0 : 
L’équation différentielle s’écrit : 


ÿ + wq y = 0 

La solution s’exprime dans ce cas : 

y(t) = A cos(cüo t + (j)) 

Cette solution est appelée solution harmonique car y(t) est une fonction sinusoïdale du temps, 
de pulsation uq , de période Tç> = 2tt/üüq et dont l’amplitude A et la phase initiale 4> son t 
déterminées par les conditions initiales : 


y 0 = A cos(0) 
yo = -w o A sin (0) 


d’où l’on tire : 


4> = — arctan 


upyp 

ÿo 



La figure suivante représente les variations au cours du temps de y(t) dans le cas particulier 
où : y 0 = 1 et ÿ 0 = 0. 
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Equations différentielles 



Oscillations non amorties 


A. 3 Equation avec second membre 

A. 3.1 Solution générale 

Soit y(t) la solution générale de l’équation différentielle avec second membre : 

ÿ + 2ôy + uly = A(t ) 

Soit unit) la solution de l’équation homogène et soit yp(t) une solution particulière de l’équa- 
tion avec second membre; r/#(t) et yp(t) sont les solutions respectives des deux équations diffé- 
rentielles suivantes : 


Vh + 2 5 ÿp + wq y h = 0 


ÿ P + 2 5 ÿ P + Wq yp = A(t) 

Les opérations de dérivation qui interviennent étant des opérations linéaires, l’addition membre 
membre des deux équations différentielles précédentes donne : 

(ÿH + ÿp) + 2 6 ( ÿ H + ÿp) + Wo {y H + y p ) = A(t) 

Ainsi la solution générale peut être obtenue en faisant la somme de la solution homogène et 
d’une solution particulière : 


y(t) = VH(t)+yp(t) 

A. 3. 2 Cas particulier où A(t) est constante 

Considérons le cas particulier important où A(t) = Aq , Aq étant une constante réelle. 
L’équation différentielle s’écrit alors : 


ÿ + 2ôÿ + ujQy — Aq 

On peut vérifier que 


constitue une solution particulière de l’équation avec second membre. Selon le cas la solution 
générale de l’équation avec second membre sera : 

- Régime fortement amorti (6 > coq) 
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y{t) = A ie s ^ + A 2 e S2t + 

- Régime critique (ô = uo) 

y(t) = (Ai + A 2 t ) e St H ^ 

u o 

- Régime pseudo-périodique (5 < wo) 

y(t) = A e~ 5t cos (unt + y) H — ° 

Dans les trois cas considérés les constantes d’intégration A\ , A 2 , A et (j) sont déterminées à 
partir des conditions initiales. Considérons à titre d’exemple le cas particulier yo = 0 et ÿo = 0 
. La résolution des systèmes d’équations obtenus permet de calculer les constantes d’intégration 
et de tracer les variations de y(t) pour les trois régimes. Les variations correspondantes de y(t) 
sont représentées par la figure ci-dessous. 



A. 3. 3 Cas particulier où A(t) = A 0 cos (Vit) : 

Nous pouvons vérifier que yp(t) = Yq cos(Dt + 9) constitue une solution particulière de 
l’équation différentielle avec second membre à condition que l’amplitude Yo et la phase 9 vérifient 
la relation : 


[cjq — D 2 ] Yo cos(f ït + 9) — 2 ô fl Yo sin (Qt + 9) = Aq cos [Vtt) 
Le développement des termes en cosinus et en sinus permet d’obtenir : 


Y 0 = 


Aq 

^[u 2 0 - n 2 ] 2 + 4 ô 2 n 2 


9 = — arctan 


OJ 

La solution complète s’écrit alors suivant le cas 
- Régime fortement amorti (5 > ujq) 


2ÔÜ 


y(t) = A\e Slt + A 2 e S2t + Yo cos (fit + 9) 
- Régime critique (ô = cuq) 
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y(t) = ( Ai + A- 2 t) e St + Yq cos (EU + 9) 

- Régime pseudo-périodique (5 < uq) 

y(t) = A e~ St cos(u;£)f + <£>) + Yq cos (fit + 9) 

Comme dans le cas précédent les constantes d’intégration A\ , A 2 , A et 9 sont déterminées 
à partir des conditions initiales. La figure suivante représente le résultat obtenu dans le cas 
particulier où : ( yo = 0 , ÿo = 0 ) . 


Régime 



Cas où A(t) est une fonction sinusoïdale 


On remarque que dans tous les cas la solution homogène tend zéro lorsque t augmente, et la 
solution générale s’identifie alors avec la solution particulière : 

y(t) a yp(t) 

Pour cette raison la solution de l’équation homogène est appelée solution transitoire tandis 
que la solution particulière est appelée solution permanente. 

A. 3. 4 Cas où A(t) est une fonction périodique du temps 

Introduction 

Nous avons étudié dans le paragraphe précédent la solution de l’équation différentielle dans le 
cas harmonique, c’est-à-dire lorsque le second membre est une fonction sinusoïdale de la variable 
t. En considérant le cas de fonctions périodiques, nous procéderons à une généralisation du cas 
harmonique. 

Développement en série de Fourier d’une fonction périodique 

Les fonctions périodiques Une fonction f(t) est dite périodique, de période T, si pour tout 
t, f(t + T) = f(t) où T est une constante positive. La plus petite valeur non nulle de T est appelé 
la période de f(t). 

Exemples : 

1. La fonction sin(7rt) reprend les mêmes valeurs pour t = 2,4,6 puisque sin[7r(f + 2)] = 
sin[7r(f + 4)] = sin[7r(f + 6)] = sin(t). Cependant 2 est la période de sin(7rf). 

2. La période de sin (Qt) ou de cos (Qt) est 2-k /EL. La période de sin (nELt) et de cos (nElt) est 
2tt / nEL. 

3. Une constante a pour période n’importe quel nombre positif. 
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4. Les figures ci-dessous représentent deux exemples de fonctions périodiques non sinusoïdales. 



Exemples de fonctions périodiques 


Définition des séries de Fourier Soit /(t) une fonction périodique de période T, c’est-dire 
telle que f(t) = f(t + T), la série de Fourier ou le développement de Fourier qui correspond /(t) 
est définie par : 

OO 

/W = f + £« n cos(nf U) + b n sin (nQt) 
n= 1 

où Q = 2ir/T, n est un entier positif, a n et b n sont les coefficients de Fourier. Ces coefficients 
sont définis par les expressions suivantes : 


1 f T 

a n = — f(t) cos (nüt)dt 

1 J o 


1 f T 

b n = — / fit) sin (nttt)dt 

1 J o 


La valeur clq/2 représente la valeur moyenne de f(t) sur la période T. Le terme sinusoïdal de 
pulsation fîi = Q , la plus faible est appelé fondamental tandis que les termes de pulsation 
Q n = nSl plus élevée sont appelés les harmoniques. 


Cas des fonctions paires et impaires 

- Une fonction est dite paire si /(— t) = f(t). Dans le cas de développement en série de 
Fourier d’une fonction paire, il n’y aura que les termes en cosinus et parfois une constante 
qui représente la valeur moyenne. Il est aisé de montrer que : 

b n = 0 




4 

T 



fit ) COS 


2irnt 

T 


dt 


- Une fonction est dite impaire si /(— t) = —fit). Seuls les termes en sinus peuvent très 
présents dans un développement en série de Fourier d’une fonction impaire. Il est également 
aisé de montrer que : 


ü n — 0 

j 4 f T / 2irnt\ 

b n = - fit) sin \-f~ J dt 
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Equations différentielles 


Solution de l’équation différentielle La fonction A(t) étant périodique, de période T, son 
développement de Fourier s’écrit : 

OO 

m = y + E» n cos (nClt) + b n sin (nQt) 

n = 1 

L’équation différentielle s’écrit alors : 


ÿ + 25ÿ + ujQy— — + E“ n cos (nü,t) + b n sin (n£lt) 


n= 1 


La réponse permanente ( ou stationnaire ) qui s’identifie avec la solution particulière, pour t 
suffisamment élevé, peut alors être calculée pour chacune des composantes de l’excitation : üq/2, 
a n cos(nflf), b n sin(rdït). On obtient alors par superposition : 


»w = ^ + E 


2Ü§ — , 

u n= 1 


&n cos + On} “h bfi 0 ^ 

V(^-"o) 2 + 4^1 
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Compléments de mathématiques 

B.l Fonctions trigonométriques 

Exercice 1 : Un déplacement harmonique est décrit par x(t) = 10cos(7rt/5) (x en mm, t en 
secondes et la phase en radians). Déterminer : 

1. la fréquence et la période du mouvement ; 

2. l’amplitude du déplacement, de la vitesse et de l’accélération ; 

3. le déplacement, la vitesse et l’accélération aux instants t = 0 s et t = 1.2 s. 

Exercice 2 : Un accéléromètre indique que l’accélération d’un dispositif mécanique est sinusoï- 
dale de fréquence 40 Hz. Si l’amplitude de l’accélération est de 100 m.s^ 2 , déterminer l’amplitude 
du déplacement et de la vitesse. 

Exercice 3 : Un mouvement harmonique est décrit par x(t ) = X cos(100f + </>). Les conditions 
initiales sont x(0) = 4.0 mm et x (0) = 1.0 m/s. 

1. Calculer X et <p. 

2. Exprimer x(t) sous la forme x = Acos(u;t) + B sm(ut) et en déduire les valeurs de A et B. 
Exercice 4 : Montrer que 

x(t) = 2 sin(u;t) + 3 cos (u;t) 

peut se mettre sous la forme :x(t) = X cos(tcf + a). Quelles sont les valeurs de A et a? 

Exercice 5 : Représenter graphiquement les variations au cours du temps d’un mouvement 
périodique décrit par : 

1. x(t) = 3sin(27rt) + 3sin(207rf) ; 

2. x(t) = 3sin(27rf) + 3sin(2.27rf). 

B. 2 Séries de Fourier 

Exercice 6 : 

1. Représenter graphiquement chacune des fonctions suivantes dont la période est égale à T. 
Préciser la parité de chacune de ces fonctions. 

(a) 

/(£) = t pour — T/2 < t < T/2 
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Compléments de mathématiques 


(b) 

J f(t) = —t pour 
\ f{t) = t pour 

- T/2 < t < 0 
0 < t < T/2 

(c) 

( f(t) = —a pour 
\ f(t) = a pour 

— T/2 < t < 0 
0 < t < T/2 

(d) 

f(t ) = t 2 pour 

- T/2 < t < T/2 

(e) 

J f{t) = 0 pour 
\ f(t) = t pour 

— T/2 < t < 0 
0 < t < T/2 


2. Calculer les coefficients du développement en série de Fourier de chacune de ces fonctions. 

3. Quelle est la valeur moyenne sur une période de chacune de ces fonctions ? 

B. 3 Nombres complexes 

Exercice 7 : Exprimer les nombres complexes suivants sous la forme z = \Z\ e t<p 

1. z = 1 — i y/3 

2. z = -2 


4. z = 5i 

5 - 2 = phf 

6. z = ( y/3 + i) (3 + 4 i) 

7 ~ \/3 — i 

'• Z ~ 3-4 i 

8. z = (2 i) 2 + 3i + 8 

B. 4 Equations différentielles 


Exercice 8 : Calculer et représenter graphiquement la solution de l’équation différentielle 

homogène : 

x + 4 x = 0 

, pour les conditions initiales suivantes : 

1. x(0) = 1 et x (0) = 0. 

2. x(0) = 0 et x(0) = 2. 

Exercice 9 : Pour les conditions initiales suivantes : 

1. x(0) = 1 et x (0) = 0. 

2. x(0) = 0 et x(0) = 2. 

3. x(0) = 1 et x(0) = 2. 

calculer et représenter graphiquement les solutions des équations différentielles homogènes sui- 
vantes : 

1. x + 5x + 4æ = 0 
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2. x + 4x + 4x = 0 

3. x + 0.3x + 4x = 0 

Exercice 10 : Pour les conditions initiales x(0) = 0 et x (0) = 0, calculer et représenter gra- 
phiquement la solution générale de chacune des équations différentielles inhomogènes suivantes : 

1. x -{- 4x = 5 

2. x + 5x + 4x = 5 

3. x + 4x + 4x = 5 

4. x + 0.3x + 4x = 5 

5. x 4x = 5 cos(3t) 

6. x + 4x = 5cos(2.1f) 

7. x + 5x + 4x = 5 cos(3 t 

8. x + 4x + 4x = 5 cos(3 1) 

9. x 4- 0.3x + 4x = 5 cos(3i) 

Exercice 11 : Calculer la solution particulière de chacune des équations différentielles inhomo- 
gènes suivantes : 

1. x + 4x = f(t ) 

2. x + 5x + 4x = /(i) 

3. x + 4x + 4x = f(t) 

où /(i) est une fonction de période T, définie par : 

f f(t ) = —a pour — T/2 < t < 0 
\ f(t ) = +a pour 0 < t < T/2 

a étant un nombre réel positif. 


H. Djelouah 



108 


Compléments de mathématiques 


H. Djelouah 



Annexe C 


Dynamique du solide en rotation 
autour d’un axe 


Exercice 1 : 

Soit un cylindre droit, plein et homogène, de rayon R, de hauteur h 
et de masse M. 

1. Trouver son moment d’inertie, par rapport à son axe de révo- 
lution (A). 

2. Utiliser le théorème d’Huyghens pour trouver le moment 
d’inertie de ce cylindre par rapport à l’une de ses génératrices 

(AO- 

Exercice 2 : 

Soit une barre homogène de faible section, de longueur L et de masse 
M. 

1. Calculer son moment d’inertie par rapport à un axe perpen- 
diculaire passant par son milieu (A). 

2. Calculer son moment d’inertie par rapport à un axe perpen- 
diculaire passant par une de ses extrémités (AO- 

Exercice 3 : 

Une barre homogène de faible section, de masse M et de longueur 
L, tourne sans frottement dans un plan vertical autour d’un axe 
horizontal (A) fixé à l’une de ses extrémités. 

1. Calculer son énergie mécanique. 

2. En déduire l’équation différentielle qui régit son mouvement. 

Exercice 4 : 

Un cylindre de masse M et de rayon R roule sans glisser sur un plan 
incliné faisant un angle 8 avec l’horizontale. 

1. Calculer son énergie mécanique. 

2. En déduire son accélération. 

Exercice 5 : 

Etablir l’équation différentielle qui régit le mouvement d’un cylindre 
de masse m et de rayon r qui roule sans glisser dans une demi 
cylindre de rayon R. On se limitera à l’étude du cas particulier du 
mouvement dans un plan vertical. 
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